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ABSTRACT 


An empirical method is developed for predicting the modal damp- 
ing of a combined parallel- stage Shuttle model by means of damping 
measurements performed on the individual substructures. Correlations 
are first determined for each component in terms of damping energy as 
a function of peak kinetic energy and modal amplitude. The results are 
then used to predict component damping energies corresponding to the 
respective kinetic energies and amplitudes that occur for the new modes 
of the combined System. Modal characteristics for the System, other 
than damping, are obtained by a real eigenvalue solution of dynamic 
equations developed by Hurty 1 s procedure of substructures. System equa- 
tions, which include component modal damping, are also solved by a 
complex eigenvalue approach for comparison with results of the empirical 
method. 


The experimental model consists of two beam-like substructures, 
each of which contains a series of discrete rigid masses, and can be 
joined into a single system through three pin-joints. Variable damping 
is provided by means of multiple paddles acting in oil reservoirs. Joint 
damping is also provided by means of a magnetic loop which is modulated 
by an active feedback circuit. Components are tested in pin- slip and 
free-free support conditions, and the System is tested in the free-free 
condition. A variety of damping and mass configurations are included. 
The empirical method is found to provide damping predictions within 
10% to 20% error, while the complex eigenvalue results deviate by as 
much as 300% ! Recommendations are given for application of the 
empirical method to a four -component, more representative Shuttle 
configuration. 
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PRINCIPAL NOTATION 
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Area 

Coordinate transformation matrix 

System coupling matrix 

Damping matrix 

Shear modulus 

Eigenvector matrix 

Damping energy dissipated per cycle 

Link length, etc. (see Figure 4) 

Modulus of elasticity 

Column vector of generalized forces 

Frequency (Hz) 

Critical damping ratio 

Moment of inertia 

Identity matrix 

Stiffness matrices 

Generalized stiffness matrix 

Element length or reference length 

Mass matrices 

Generalized mass matrix 

Denotes normal mode elastic quantities 

Poisson's ratio 

Column vector of nodal forces 

Column vector of generalized coordinates 

Column vector of generalized applied forces 

Column vector of displacements or generalized 
coordinates 

Denotes radius or rigid-body quantities 

Kinetic energy 

Tube thickness or time 
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PRINCIPAL NOTATION (Cont'd. ) 


Internal energy 

Column vector of nodal displacements 
Modal displacement matrices 

Peak modal amplitude for steady-state sinusoidal 
vibration 

Longitudinal, lateral, and rotational displacements, 
respectively, at the node points. 

Phase angle 

Ratio of forcing frequency to j-th natural frequency 
Circular frequency 

Denotes quantities are with respect to coordinate 
system where all nodes are interior nodes. Also 
denotes normalized quantity. 

Denotes quantities are with respect to coordinate 
system where some nodes are on boundary and the 
rest are interior nodes 
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I. INTRODUCTION 


Prediction of damping in complex structures is recognized to be 
extremely difficult because of the presence of various damping mechanisms 
and the current lack of knowledge in synthesizing even the most elementary 
mathematical models into combined structural systems. As a result, in 
initial design procedures in the past, damping characteristics of aerospace 
systems have usually been predicted at the modal level by empirical methods 
based on data derived from previous vehicles^)*, and when possible, ulti- 
mately verified by actual measurements of modal damping in full-scale 
dynamic tests. Needless to say, such full-scale testing of large systems is 
extremely difficult and expensive at best, and has been deemed unfeasible 
for the Space Shuttle, which will be the largest and most complex aerospace 
structural system yet attempted. 

In order to provide information necessary for predicting dynamic 
response for Space Shuttle, it is recognized that some new procedure must 
be developed for this multi -component system. As a reasonable approach, 
full-scale testing of the Booster, Orbiter, and any other major components 
appeared possible, but a requirement then exists for synthesizing this in- 
formation into a prediction of response for the combined System. Substruc- 
turing methods are available^), and are presently being further developed^) 
for handling prediction of dynamic characteristics of undamped structures. 
However, no accurate method has been derived for including damping in 
such structural synthesis procedures. Therefore, the purpose of this report 
is to present the development of a method which shows great promise for 
accomplishing this objective. The philosophy followed in the research pro- 
gram leading to these results has been to perform dynamic tests on individual 
substructures of a two- component, parallel-stage, rather simple Space 
Shuttle model, develop a mathematical model for the components, combine 
all of the information into a mathematical model for the combined System, 
and compare predicted dynamic characteristics with those obtained from 
further tests of the combined System. 

At the outset of this program it was recognized that enormous ef- 
forts have been exerted in the past^’ ^ to develop mathematical models of 
various elementary damping mechanisms, and this work will continue be- 
cause of the great uncertainty present in the current state-of-the-art in 
combining these mechanisms in complex structures. It was also recognized 


Superscript numbers in parentheses denote references on page 92 
of this report. 
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that the most reliable methods currently available for measurment of 
typically light damping in complex structures having reasonable modal 
separation, are those which determine an overall modal damping at 
resonance or free decay from resonance conditions. The beauty of the 
concept of modal damping, of course, is that it results in a single 
number that represents a reasonable approximation to the aggregate 
contributions of a multitude of complex damping mechanisms internal 
to the overall structure. Once it is determined, and some form of relatively 
simple equivalent mathematical model is matched with the modal values 
of mass, stiffness, and damping, prediction of dynamic response for that 
component will be satisfactory, at all intermediate frequencies, since 
damping forces dominate only in the vicinity of resonances. Nonlinearity 
of damping with amplitude, and inevitable relatively wide scatter in 
damping measurements are two unfortunate dilemmas which must often 
be faced with this otherwise reasonable approach to the dynamic predic- 
tion problem. 

In view of the success of modal concepts for synthesis in the past, 
it appeared reasonable to consider their use for Space Shuttle predictions 
also, although it is recognized that considerable complications result 
because of the combination of the substructures into an even more complex 
System. System natural modes can be expected to occur at frequencies 
which are intermediate to those for the individual components, and as a 
result, damping forces which were insignificant in the dynamic response 
at the intermediate frequencies for the individual structures, now become 
dominant for the new modes and frequencies. Therefore, the prediction 
method outlined herein, is based on determination of modal damping 
energy for each component from its own resonance tests, extrapolation 
of the contribution of each respective component to the total modal damp- 
ing energy at the new frequencies, and summation of the results for all 
components. Results for the method will be shown to be far superior to 
others that were considered. Limitations will be discussed, along with 
additional work recommended for further verification. We begin by pro- 
viding a description of the physical model used in the study, proceed with 
derivation of mathematical models of the System, continue with a detailed 
development of the damping synthesis method, and finally compare results 
for a multitude of tests and predictions for modal damping and dynamic 
response of a variety of damping and mass configurations. 
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II. DESCRIPTION OF PHYSICAL, MODEL 


A. Configurations 

A conceptual drawing of the Space Shuttle model is shown in 
Figure 1. It consists basically of Booster and Orbiter model joined 
together with a pin-joint near the top, and two pin-joints at the end of 
a rigid link near the bottom. The lower joint simulates a slip-joint for 
all practical purposes. Each component was made of a 2024 aluminum 
tube of 1-inch O. D. and 0. 035 -inch wall. Multiple rigid masses were 
positioned on each to represent the two different mass configurations 
indicated in Figure 1. Approximate locations of the masses can be 
scaled from the figure. These masses were made of split disks of 1/2-inch 
thick mild steel and all were nominally 6 -inch diameter and weighed 
approximately 3-3/4-lb each. Those masses attached to connecting joints 
included cutouts where necessary to accommodate the joints. These 
pin- joints were provided through the use of 1 /4-inch diameter, Bendix 
Type 5008-600 flexure pivot bearings, that are designed to have ex- 
tremely low inherent damping. 

A list of the numerous configurations to be considered is given 
in Table 1(a), while an explanation of the notation is given in Table 1(b). 

It can be seen that two mass configurations were included for several 
different damping configurations which will be explained momentarily. 

Both pin- slip and free-free boundary conditions were investigated for 
the components in order to determine what effect this would have on pre- 
dictions for the System. (This consideration would also help indicate 
the type of support structures required for full-scale component testing. ) 
Pin- slip component tests were conducted by mounting the substructure 
directly to a rigid wall through the upper pin and lower slip- joint. Free- 
free component tests were conducted by hanging the substructure from a 
low frequency suspension which was attached to eye -bolts at each side of 
the second-from-the-top mass (this point was estimated to be a node for 
most of the lateral bending modes). Finally, System tests were conducted 
by having the two components attached together with a pin-joint near the 
top and simulated slip-joint near the bottom, while the entire System was 
supported on a low frequency suspension attached to a rigid knife-edge 
nearly under the upper pin-joint. Thus, a simulated overall free-free 
suspension was employed. A photograph of the model in this configura- 
tion is shown in Figure 2. 





Pin -Slip Components 


1. B-P-L1-D0-M1 

6. O-P-Ll-DO-Ml 


2. B-P-L1-D1-M1 

7. O-P-Ll-Dl-Ml 


3. B-P-L1-D2-M1 

8. O-P-Ll -D2-M1 


4. B-P-L1-D3-M1 

9. O-P-Ll -D3-M1 


5. B-P-L1-D2-M2 

10. 0-P-L1-D2-M2 


Free-Free Components 



11. O-F-LO-DO-Ml 

16. B-F -L0-D0-M1 


12. O-F -L0-D2-M1 (RB Incl.) 

17. B-F-L0-D2 - Ml 

(RB Incl.) 

13. O-F -L0-D3-M1 (RB Incl.) 

18. B-F -L0-D3-M1 

(RB Incl.) 

14. O-F-LO-Dl - Ml (RB Only) 

19. B-F -L0-D1 - Ml 

(RB Only) 

15. O-F-L0-D2-M2 (RB Only) 

20. B-F-L0-D2-M2 

(RB Only) 

Additional Connecting Link Damping 

Included 


21. LI 

23. L3 


22. L2 

24. L4 



Combined System 


25. 

S-F-L1-D0-M1 

30. S-F -LI -D2/D3-M1 

26. 

S-F- LI -D2-M1 

31. S-F- L2-D2/D3-M1 

27. 

S-F -L2-D2-M1 

32. S-F -L3-D2/D3-M1 

28. 

S-F -L3- D2-M1 

33. S-F-L4-D2/D3-M1 


29. S-F -L4-D2-M1 

34. S-F-L1-D3-M1 

35. S-F -LI -D1 -Ml 

36. S-F-L1-D2-M2 




TABLE 1(b) 
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CONFIGURATION NOMENCLATURE 


Vehicle Identification 


B - 

Booster 


O - 

Orbiter 


S - 

Combined System 


Support 



P - 

Upper Support Pinned, Lower Support Slip 
or Double Pin 

F - 

Components are Free-Free, System ha 
Free -Free Support but components are 
joined together as with P- Support 

,s 

Connecting 

Links 


LO 

- No Connection 


LI 

Very Low Damping (Three Flexure Pins) 

L2 

- 0. 5% Critical Damping for System in 

First Mode 

Air 

L3 

1. 0% Critical Damping for System in 
First Mode 

Air 

L4 

2. 0% Critical Damping for System in 
First Mode 

Air 

Overall Structural Damping (Fluid Average Viscosities in cp) 


Booster 

Orbiter 

DO 

In Air 

In Air 

D1 

230 

600 

D2 

1070 

1080 

D3 

4330 

4900 

Mass Distribution 



Ml | 


M2 


See Figure 1 






B. 


Damping Mechanisms 


It was desirable to use several levels of damping for the overall 
structure, as well as provide joint damping in the System. Since re- 
peated tests were contemplated, it was necessary to provide damping in 
a form that could be controlled as accurately as possible, and be as 
repeatable and reliable as possible. The actual mechanism whereby 
damping was provided was arbitrary, so long as it had the above charac- 
teristics, and preferably provided damping forces which were linear with 
amplitude for the sake of increasing measurement reliability. After 
consideration of several mechanisms for structural damping, a series 
of wipers acting in oil reservoirs was selected, and the solid structure, 
including pin-joints, was designed to have as low inherent damping as 
possible. 

Angle-shaped wipers were attached to the masses indicated in 
Figure 1. Oil reservoirs supported in their gantry-like structure can 
be seen in Figure 2. Thus, the model itself was suspended from the 
laboratory ceiling, and the dampers touched only the oil in its space - 
fixed reservoirs. Oil viscosity was monitored along with its tempera- 
ture to allow small corrections to damping values where appropriate, 
as will be explained later. Average values of these viscosities are 
shown in Table 1(b). The D1 and D2 damping was provided by standard 
motor oils, while the D3 damping fluid consisted of a mixture of heavy 
motor oil and STP oil additive. 

Connecting link damping was provided on the Booster side of the 
lower slip-joint according to the approximate values given in Table 1(b). 
Figure 3 shows the mechanism which was found to be most reliable. It 
consists of a permanent magnet whose flux path is split by a small gap. 
Part of the magnet is attached to the rigid link while the other part is 
attached to the connecting rigid mass of the Booster. Relative motion 
between the two parts of the System is measured by accelerometers A x 
and Aj, whose output is converted to relative velocity and fed back as a 
voltage into a coil which proportionately modulated the flux path. The 
low pass filter was necessary to avoid high frequency instabilities in the 
feedback loop. Although only one such damper was used, similar damper 
could have been employed at the other two pin-joints, and the methods to 
be developed and results obtained would still be valid. 
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Figure 3. Link Damping Mechanism 
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III. DEVELOPMENT OF ANALYTICAL MODELS 


A. Coordinate System 

The coordinate system used to describe the motion of the compo- 
nent Booster and component Orbiter acting as individual substructures, 
as well as the motion of the combined System, is shown in Figure 4. 
Specific values for masses, lengths, etc. , are given for various con- 
figurations in Appendix A. At the i-th interior node point, which is 
located at the center -of- mass of the i-th concentrated mass, three 
planar degrees -of-freedom (x t , z t , 9 1 ) are permitted with the positive 
directions taken as indicated. Rigid body coordinates are referred to 
a node located at the center-of-mass of the component or System. In 
addition to these nodes we define a set of boundary nodes denoted by PI 
and P2 at the pin-joint interface and by SI and S2 at the slip-joint 
interface, respectively. Three planar degrees -of-freedom are also per- 
mitted at these nodes, until boundary constraints are applied later. 
Selection of boundary nodes in this manner facilitates the solution of the 
component eigenvalue problem since the boundary conditions may be speci- 
fied at these nodes. Likewise, when deriving the transformation matrices 
which couple the components together it is again convenient to specify 
displacement interface continuity in the appropriate directions at these 
nodes. 


Assuming zero compliance between, say, nodes 2 and PI on 
the Orbiter, the relationship between the boundary displacements 
(node PI) and the interior displacements (node 2) is given by 


X l,l = X 2 + 6 1 9 2 


or in matrix form 


f v 

x l,l 


1 

0 



x 

-*■ 2 

< z l.l 

> = 

0 

1 

0 

< 




0 

0 

1 


CD 

(0 


( 1 ) 

( 2 ) 

(3) 


(4) 


Similar relationships apply at the node pairs (5, SI), (8, P2), and (11, S2), 
respectively. The additional coordinates introduced at the boundary nodes 
are not independent but are related to the interior coordinates by relation- 
ships of the type given by Equation (4). Therefore, for free-free 
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components, the total number of Booster degrees-of-freedom is 7X 3 = 21 
and the number of Orbiter degrees-of-freedom is 6X3= 18. 

In order to formulate the equations of motion for the substructures 
it is necessary to establish some basic relationships which will be used 
throughout the remainder of this report unless otherwise noted. Let each 
triplet of independent coordinates (x 1( z 6j) corresponding to the in- 
terior nodes be denoted by a matrix vector of the form 



( 5 ) 


For either the Booster or the Orbiter the complete displacement vector 
may then be written as 


' u. 


(u3 = i 


U : 


■ 


u 


m / 


( 6 ) 


where m = 7 for the Booster and m = 6 for the Orbiter. 


In general, the total displacement at a node will be the sum of 
three displacements: 

(a) A displacement component corresponding to a rigid-body 
rotation and rigid-body translations 

(b) A displacement component corresponding to constraint 
modes in the case where redundant constraints are present 

(c) A displacement component corresponding to the natural 
or normal elastic modes of the structure. 

Since the constraints for our Booster and Orbiter models are statically 
determinant for both the pin- slip and free-free cases, modes associated 
with redundant constraints are not applicable. Therefore, the total dis- 
placement at a node is comprised solely of rigid-body and/or elastic 
normal mode displacements. In light of this. Equations (5) and (6) may 
be written as 





( 7 ) 
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and 


(*xl 


D 


- N 
U 2 

* > 
• 

U 2 
< ; 

>■ I 

U 2 

• 

• 

\ 





• 

- N 
U 

V B 


or 

U3 = {u R } + {u N } 


where 


( 8 ) 


(9) 


R denotes rigid-body displacement 

N denotes normal elastic mode displacement 

The displacements on the right-hand side of Equation (9) may be 
defined in terms of a set of normalized displacement functions or modes 
and a set of generalized displacements. Thus 

{u} = [0 R ] [p R 3 (10a) 

{u N 3 = [0 N ] Cp N 3 (10b) 

where the element 0 tJ in the modal matrices is the displacement at 
point i in the j-th mode. If these equations are substituted into 
Equation (9), the total displacement column may be written as 

[u3 = [0 R ] {p R 3 + [0 N ] {p N } = [ 0 ] {p} (ii) 


where 

[0] = [0 R ;0 N ] 


( 12 ) 


and 


{p} = 



(13) 


From Equation (4) we note that for our model the boundary dis- 
placements may be expressed in terms of the interior displacements at 
a given node. Conversely, if we select the boundary displacements as 
the independent variables then we may write 
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or 



{u 2 } - [ a lfl ] {u x>1 } 


(14) 


(15) 


At the other boundary nodes we have 

CM = [ a 3,x] 

^8 J = [“1,2] 

^■ U ll^ = [ a 2,s] 


similarly the relationships 


^2,X 5 

(16a) 


(16b) 

^a, 2 } 

(16c) 


where [ a 2 ,xJ* . . . [a 2/3 ] may be determined from the kinematic constraints 
between the nodes. Using Equations (15) and (l6a,b, c) and the node-point 
designations in Figure 4, we may write the following relationship for the 
Orbiter 


*x 


1 

— 





“i.i 

0 


^1,1 



1 







1 



u 5 


0 

tt 3,X 


^3,1 

1 

to 
13 
1 

0 


I 

0 

j*6 


(17) 


or similar to Equation (9) 

£u} 0 = [H £u} 0 = [a] 0 ({u*} 0 + {u N } 0 ) 
Also for the Booster we have 


— — 


— — 


— " 

u 7 


I 


u 7 

^8 


°i,a 0 


U X,2 

31 

CO 


I 


*9 

S 10 

= 

I 


5 10 

S X1 


a 3,3 


^3,3 

*13 


0 I 


*X3 






__ 

U i3 

B 

I 

B 

U. — 


(18) 


( 19 ) 


1=3X3 identity matrix 
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or again similar to Equation (9) 

tu} B = [a]g {u } 8 = [a^ ({u r } b + {u N } B J (20) 

The matrices {u} 0 and {u} 0 now include displacements at the two 
boundary nodes as well as at the n-2 interior nodes (n = number of 
interior nodes for the component). 

Equations (18) and (20) can be expressed as one equation with the 
understanding that subscript B applies to the Booster only and that sub- 
script O applies to the Orbiter only. Hence, 

^B,0 = r a i,0 < 21 ) 


For both the Booster and Orbiter these displacements can also be ex- 
pressed in terms of a set of normalized displacement functions and a set 
of generalized displacements. Thus 


{u R l = [0 R ] {p R 3 (22) 

OT = [0 N ] Cp N 3 (23) 


For either the Booster or Orbiter we then have from Equation (20): 



in,, (rA.„ {A,. + [*"i,. <A,.) 

= f a i,o ^p4,o 


(24) 


where again as in Equations (12) and (13) 



(25) 


and 



( 26 ) 


For the free-free eigenvalue problem it is simpler to use the {u} coordi- 
nate system and for the pin- slip eigenvalue problem the {u} coordinate 
system. A development of the component eigenvalue problems follows. 
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B. Component Eigenvalue Problems 


1. Pin- Slip Configurations . A finite element approach was used 
in deriving the stiffness matrix for each component. The component was 
discretized into 6 elements for the Booster and 5 elements for the Orbiter 
with the length of each element taken as the distance between the interior 
nodes as shown in Figure 4. A derivation of the element stiffness matrix 
for an arbitrary element with length is given in Appendix B. 

Assembly of the element stiffness matrices results in the global stiffness 
matrix [K] for the component. The concentrated masses and moments 
of inertia as well as the distributed mass of the tube are lumped at the 
interior nodes. This results in a global mass matrix [M] for the compo- 
nent. Since these stiffness and mass matrices are expressed solely in 
terms of the interior coordinate system we now include the boundary nodes 
by equating the internal and kinetic energies of the components in the two 
systems [Equations (11) and (24)]. For either component and mass con- 
figuration, equating strain energies in the two systems gives 

^ B ,0 = 2 |_ uJb »0 ^^ 8,0 = 2 L 5 Jb ,0 [^ 1 b ,0 ^^,0 (^ 7 ) 


Substituting Equation (21) into the right-hand side of Equation (27) gives 

2 lAk ° 0 ^ U ^ 8,0 = 2 [_ u ] B ,0 f a 3 s » 0 0 [ a L ,0 ^ U \,0 (^ 8 ) 

from which we obtain the relationship 

[*1,0 = K,o [*l,o [*1,0 ( 29 ) 


Similarly by equating the kinetic energies in the two systems we have 

[“1,0 = K, 0 [al >0 (30) 

For both the Booster and Orbiter component pin- slip configurations 
the boundary conditions at the joints (boundary nodes) are as follows: 


Booster: 

X 1 >2 - 2 1,2 ~ 0 

at pin -joint 

(31a) 


2 a, 2 = 0 

at slip-joint 

(31b) 

Orbiter: 

x l,i = 2 1,1 = 0 

at pin-joint 

(32a) 


2 2,1 = 0 

at slip-joint 

(32b) 


The stiffness and mass matrices are now partitioned into the following 
submatrices: 
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t K Lc 


I M L< 


K„ , K 


fc 


*:■ t 


L K cf ! 

KccJ 

* A i 

M ,, 

M, 

ff 1 
1 

t C 

A 1 

~ 7 ~ 

L M cri 

M 

CO 


B >0 


(33a) 


(33b) 


The elements in the displacement matrix (u} B>0 are also partitioned so 
that the elements corresponding to the constrained boundary displace- 
ments [Equations (31a, b) and (32a, b)] form one submatrix and the free 
displacements form another submatrix. Thus 



(34) 


where 

f denotes the free displacements 

c denotes the constrained displacements 

— > means replaced by. 


The equations of motion for the undamped component undergoing free 
vibration are therefore given by 



(35) 


Noting now that the submatrices { q c ]^ , 0 and {q 0 ^, 0 are zero 
due to the boundary conditions, it follows from Equation (35) that 


~ ^^b»o (36a) 

[^4 0 ft»o + [Kcfi.o = f^o^B»o (36b) 


Equation (36b) gives the boundary reaction forces once the displacements 
at the free coordinates are known. Assuming simple harmonic motion of 
the form 



18 


U f i,o 


Cq°} B ,c 


itut 

e 


(37) 


and substituting into Equation (36a) we obtain the standard eigenvalue 
problem for the pin- slip boundary conditions 


[I] - P] 


= 0 


(38) 


where 


[D] = [M [K, t ]„ 0 


(39) 


A modified Jacobi method was used to obtain the eigenvalues and mode 
shapes corresponding to Equation (38). The mode shapes are normalized 
to unity displacement at the maximum displacement in each mode and are 
coupled in the axial and transverse directions due to the pin- slip boundary 
conditions. Only the transverse displacements are given in this report, 
a detailed description of which is given in Section V. Since all rigid-body 
degrees-of-freedom have been suppressed, the modes obtained are the 
elastic natural modes of the component. 


The generalized stiffness and mass matrices may now be written as 


[*f NN ] = [0 N f ; 0 " = o] 


n : 


B,0 


where 


[#?] 

0 


8,0 


A | A 


!>, 1 

K ff .K, c 


A 1 A 


00 = 0 
B,' 

|K c ,|K e J 

B, 

0 


(40) 


is a matrix of normalized modal displacements for 
the unconstrained coordinates 

is a null matrix since the displacements are zero 
at the constrained coordinates. 


and 


[^ NN ] =[0 N f ;0 N o = O] 

B,0 B,0 


i 


M 0 f'M c0 . 


f; 


B.O 


L0= = O. 

B ,0 


(41) 


2. Free- Free Configurations . The component free-free con- 

figuration eigenvalue problem can be set up in two different ways which 
we denote as the direct method and the transformation method. Both 
methods are required for subsequent combined System analysis, and 
will be included here to provide an indication of numerical accuracy in 
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component analysis. The transformation method requires a knowledge 
of the pin- slip modes in order to determine the free-free modes. Since 
in the pin- slip configurations the rigid-body modes have been suppressed, 
we must now synthesize these modes into the equations of motion. This 
is necessary for both the component case and for the System case where 
pin- slip modes are used as generalized displacements. 

a. Direct Method. In the direct method the stiffness and 
mass matrices [k] and [m] are obtained in the same manner as out- 
lined in Section B. 1 for the pin- slip configurations. These matrices are 
written in the {u} coordinate system with the displacements defined at 
the component interior nodes. Since there are no boundary constraints 
for the free-free configurations the component undamped equations of 
free motion are 


[M] b>0 Cu 3 b , 0 + [K] b>0 fu} B , 0 = {0} 


(42) 


Assuming simple harmonic motion of the form 

iOJt 


tulb.o = £u°k, 0 e' 


(43) 


and substituting into Equation (42) we obtain the standard eigenvalue pro- 
blem for the free-free boundary conditions 


ou 2 [I]- [D] 


= 0 


(44) 


where 


[D] = [MJ^ 0 [K] b , 0 


(45) 


A modified Jacobi method was used to obtain the eigenvalues and mode 
shapes corresponding to Equation (45). The mode shapes are normalized 
to unity displacement at the maximum displacement in each mode and are 
uncoupled in the axial and transverse directions. Since we have con- 
sidered only planar motion in this study, we obtain three rigid-body modes 
in addition to the elastic natural modes of the component. A detailed des- 
cription of the free -free component natural modes and frequencies is given 
in Section V. 


The generalized stiffness and mass matrices for free-free compo- 
nents may now be written as 



o 1 

o 


1 

DC 

1 

— 71" 


— 

Lo fr NN 


0 N 

B, 

0 

L B 5 


[k] [0 R ; 0 N ] 

B,0 B ,0 


( 46 ) 


and 
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m 


8,0 


^ RR ; o 


1 

or 

I'a 

1 



= 

-H 

or 1 
QC 

ST 

o 

i 


L 0 N J 

B, 

0 

B, 


. — a ■ — D 


[M] [ 0 , 0 J 

B ,0 B,0 


where 


(47) 


[0 R ] is the matrix of normalized rigid-body displacements 

[0 N ] is the matrix of normalized elastic modal displacements 


and both are synonomous with the modal functions defined in Equations 
(10a, b). It should be noted that all submatrices in the generalized stiff- 
ness matrix associated with rigid-body displacements are null since no 
internal work is done during a rigid-body displacement. 


b. Transformation Method . In the transformation method 
we attempt to obtain the free-free component eigenvalues and mode shapes 
by using our analytically derived pin- slip modes and by appropriately in- 
cluding rigid-body motion in the equations of motion. Recall that in 
Section IIL B. 1 we obtained the pin- slip modal matrix 


I 




(48) 


where mode j is specified in column j of this modal matrix and includes 
displacements at the boundary nodes as well as at the interior nodes. We 
now have to include rigid-body motion with the pin- slip modes. Let us 
therefore define in the {u} system three rigid-body modes in the axial-a, 
lateral-4, and rotational-r, directions, respectively, such that the rigid- 
body displacements are related to the modes and generalized displacements 
by the following relationships: 


{u a 3 = C 0 g 3 p a = 


1 

0 

_ 0 _ 

T 

0 

0 


/Pa 


i 

V o ' 


(49a) 
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{u/} = {0i} p £ = 


0 

1 

0 

o“ 

1 




0 

1 

V 0 > 


Cu r } = {0 r } p r = 


{ 0 ^ 
®r 

l n 

~ 6 " 



p 


r 


o 

•A. 

V 

V i lii 


(49b) 


(49c) 


where p a , p i, and p r are the respective generalized coordinates. In 
Equation (49c) a normalized rotational mode is selected such that a 
rotation of 9 = 1/4 about the component center-of-gravity produces a 
unit displacement at the tip of the component and 4 is the distance be- 
tween the eg and the tip node. Thus, z* is the displacement in the 
lateral direction at interior node j due to a rotation about the compo- 
nent eg equal to 9 = 1/1. For both the Booster and Orbiter, z * was 
equal to 1 at the bottom node. Equations (49a, b, c) may be combined 
so that 


{u r } = [0 a 0ji 0 P ] < px > = iy ] cp R 3 

p. 


From Equation (21) we had that 


^^b»o - [ a ]B»o ^ b , o 

If we multiply both sides of Equation (50) by [a]^ we have 

[a]B?o U R ) = [«].* y R ] Cp R 3 


(50) 


(21) 


(51) 
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or 


{u R 3 = [0 R ] {p R } 


(52) 


which conforms with Equation (22). Equation (52) now gives the rigid-body 
displacements in the {u} coordinate system. Similar to the normal mode 
matrix [0 N ] the elements in [0 R ] are arranged as indicated by Equation 
(48) so that 




(53) 


where again 


means replaced by. 


Combining Equations (48) and (5 3) into one equation we have for 
both the Booster and the Orbiter 


[0] 


[* R ; n. 


B,0 


and 


(54a) 


["]b,o = C^] B ,o Wb,o ( 54b ) 

Equation (54b) now gives the displacement at a point on the component as 
a linear combination of pin- slip modes and rigid-body displacements. The 
kinetic energy in the free-free configuration may be written as 


t b,o = i lAl 8,0 [M ] b>0 Cu } S)0 

= i LpJs'O [M] b ,o [0] b ,o 


(55) 


from which we define the mass matrix 


[^]b,o 



1 * RN 1 

m 

1 m 

A NR 

m 

1 » 

1 z 
i 


(56) 


Note here that the matrix [?7! NN ] is already available from the pin- slip 
analysis and is given by Equation (41). The matrix [m RR ] represents 
the rigid-body masses while the matrices [m RN ] and [m NR ] represent 
the coupling masses. Using the stiffness matrix derived in Equation (40) 
we may now write the equations of motion for the undamped free-free 
components as 
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~ * RR 
m 

*C,Rn" 

m 

i 

fp"| 


0 

! °" 


( p " : 


) 

* NR 

m 

i NN _ 

B.O 

|7j 

> + 
B,0 

0 


B >0 

KJ 

■ { 

B >0 

l) 


[mjg.o {p} B ,o + [£] b ,o fp} B ,o = (57) 

Assuming simple harmonic motion the eigenvalues and mode shapes 
corresponding to Equation (57) were determined by a modified Jacobi 
method for the Booster component using the first six (6) pin- slip modes 
and three (3) rigid-body modes as generalized displacements. Eigen- 
value results of the analysis are presented in Table II and compared 
with the results obtained by the direct method. 


TABLE II. COMPARISON OF NATURAL FREQUENCIES 
OBTAINED BY DIRECT AND TRANSFORMATION METHODS 


Mode 

No. 

Natural Frequency 
Direct Method 

Natural Frequency 
Transformation Method 

%A* 

1 

8.53 

8.54 

+ 0. 12 

2 

24. 08 

24. 13 

+ 0. 21 

3 

45. 72 

45. 60 

- 0. 26 

4 

75.25 

75. 32 

+ 0.09 

5 

106. 12 

106. 29 

+ 0. 16 

6 

203. 39 

— 

— 


Based on direct method. 


Eigenvector results can also be used to define the following 
generalized stiffness and mass matrices 


and 


m = 

B,0 


cp R 

cp\ 

B,0 


[k] [£ R i<P N ] 

B ,0 B,0 


m 

B,0 



B,0 


[m] [<P R |<P N ] 

B ,0 B ,0 


(58a) 


(58b) 
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where 


A p 

[ Cp ] is the matrix of normalized rigid-body displacements 

r A N 1 

[cp ] is the matrix of normalized elastic modal displacements. 

Now since the same component is represented in both the direct and 
transformation methods of analysis we must have 


Mb,0 = 

m B ,o 

(59a) 

a) 

o 

II 

m B , o 

(59b) 


Use of these equations allows a check on the numerical accuracy in the 
problem, and is particularly useful for checking whether a sufficient 
number of elastic modes have been used in the analysis. 

As can be seen from Table II, the natural frequencies obtained 
by the transformation method are in excellent agreement with those ob- 
tained by the direct method through the fifth mode. Inspection of the 
mode shapes also showed excellent agreement although they are not 
presented here. Inclusion of rigid-body modes in the pin- slip compo- 
nent mode synthesis for the combined System follows the same procedure 
as outlined in the transformation method. The results obtained from this 
example provided a good qualitative and quantitative indication of results 
to be expected in the System analysis which is discussed in Section III. D. 

C. Component Forced Response 

The component forced response analysis was based on the well- 
known assumption of modal damping properties which is valid for systems 
having light damping and well separated modes. That is, at resonance, 
the response is due principally to a single mode. More will be said on 
this matter in a later section. For the present discussion, however, we 
take the damping to be viscous and proportional to the linear critical 
damping for each component mode. From damping measurements made 
in the laboratory on each component we obtain a diagonal matrix of damp- 
ing coefficients which represent the percentage of critical damping for 
each elastic mode. In addition to the modal damping we also measure the 
damping obtained during rigid-body motion at very low frequencies, since 
in the free -free configurations this damping will be present due to the 
damping mechanism employed in the experiments. Note that for a real 
vehicle this rigid-body damping might correspond to low frequency aero- 
dynamic damping. 

Assuming linear superposition of the modal displacements in the 

form 
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fu} B , 0 = [^ N ]b,o^P N 3 s ,o pin-slip case 


£ 5} b , 0 = [ 0 R ! 0 N ] 


B >0 


free-free case 


B >0 


(60) 

( 61 ) 


we may, by orthogonality of the normal modes, write the equations of 
motion for forced response as 



which holds for both the pin- slip and free-free cases so long as the 
appropriate matrix elements are utilized. Further, all the square 
matrices on the left side of Equation (62) are diagonal and the damping 
matrix is given by 


fc- RR 

°i 


^rr 

ol 


'u) RR ' 

ol 


>“1 o’ 

~ 

CM 

II 



------ 


_ 

1 

L 0 


B ,0 


n m 

B >0 

L 0 i 

U) NN _ 

8,0 

_ o 1 Q J 


(63) 


where 


0 ) 


RR 


V n 


B ,0 


square diagonal matrix of rigid-body and 
elastic modal natural frequencies (note that 
[u) RR ] is null). 


(64a) 



l 

•NN 


B >0 


square diagonal matrix of rigid-body and 
elastic modal critical damping ratios 


(64b) 



vector matrix of applied generalized 
forces 


(64c) 


Equation (62) represents a set of m + 3 (m=number of elastic modes) 
uncoupled equations in the generalized p coordinates which is applicable 
to both the free-free and^pin-slip component configurations provided that 
the submatrices #? RR , C- R , and & are set equal to zero for the pin-slip 
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configurations. The j-th uncoupled equation of motion in the generalized 
coordinate corresponding to the j-th elastic mode may be written as 


p/ + 2 C/" U)/" p/ + to/" 2 p/ = y/ sin itA/" 


(65) 


For steady- state excitation a particular solution to Equation (65) is given 
by 


p/ = A/ sin (to t - Y/) (66) 

where to is the forcing frequency. If we define fij = to/lo/ N and substi- 
tute Equation (66) into Equation (65) we obtain for A/ and Y," 


?//(%/»* U)"" 8 ) 

V = [ ( 1 -Qf ) 2 + ( 2 C/" «,) a ]^ 


(67) 


= arc tan 


\ 1 - fi , 8 / 


( 68 ) 


The j-th uncoupled equation of motion in the generalized coordinate 
corresponding to the j-th rigid-body mode may be written as 


P 


R 

i 


+ 2 ^j 


RR 


top/ 


= y/ sin (OtA J RR 


(69) 


For steady- state excitation a particular solution to Equation (69) is 
given by 


p/ = A/ sin (tot - Y/) 


(70) 


Substituting Equation (70) into Equation (69) we obtain for A/ and Y/ 

?//(»?/ R (5 s ) 


[i + 4C/ rS 3 2 


(71) 


Y/ = arc tan (-2 Cj RR ) 


(72) 


Once the generalized coordinates have been obtained substitution into 
either Equation (60) or Equation (61) gives the desired displacements at 
the node points. Results of the component forced response analysis are 
presented in Section V for each configuration. 
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D. System Equations 

The procedure used to develop the System modes from the Booster 
and Orbiter component modes essentially follows that of Hurty(^). The 
component modes, as considered by Hurty, fall into three categories: 
rigid-body modes in which the component is displaced without deformation, 
modes which exist only if the system of constraints on the component is 
indeterminate (called "constraint modes"), and normal modes that define 
displacements relative to the constraint system. 

hi the first category, three rigid-body modes are considered for 
each component which consist of two rigid-body translations and a small 
rotation. Since the system of constraints on our model is determinate, 
the modes in the second category do not exist. In the third category we 
have normal modes that define displacements relative to a constrained 
system (pin- slip modes) and normal modes that are obtained with no con- 
straints present (free-free modes). 


1. Component Matrices . Since the procedure used to develop 
the System analysis is identical for both the pin- slip mode synthesis and 
the free-free mode synthesis, we will not repeat th& procedure separately 
for each case in this report. We would, however, like to make the 
following distinction between the two cases which should be observed in the 
development of the System equations: (1) For System synthesis utilizing 

pin- slip modes the component mass matrix is given by 


[ m] B , 0 


-RR 

m 

*rn 
i m 

i 

-nr 

m 

i i 


J B >0 


(73) 


ADD A Dll A MO m 

where the submatrices m , m , and m are given by 

[m RR ] e , 0 = [0 R ] b T , o [M] s ,o [0 R ] 6 ,o (74a) 

[m RN ] s ,o = [M] b , 0 [0 N ] b ,o (74b) 


and 


[m NR ] B , 0 = [0X O [M] b , 0 [0 R ] 


B ,0 


(74c) 


These submatrices were obtained by using the transformation procedure 
outlined in Section III. B. 2.b. Note again that the submatrix [#J NN ] is 
already available from the pin- slip modal analysis. (2) For System syn- 
thesis utilizing free-free modes the component generalized mass matrix 
is given by 
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m B, 0 


7 /\™ i 0 

i 

o 1 7k m 


8,0 


(75) 


where both the submatrices [ 77? RR ] and [% NN ] are available from the 
free-free modal analysis [Equation (47)]. With this distinction in mind, 
we proceed to develop the System equations. 


For the Booster we may write the equations of motion with equiva- 
lent viscous damping as 

[M] 0 tp} 8 + [£]„ tp) B +[*]> lp} 0 = £»<«), (76) 

and likewise for the Orb iter 

[»!]„ tp3 0 + [<J] 0 £p}„ + [)f] 0 tp) 0 = t»(t) }„ (77) 


where for pin- slip mode synthesis the mass matrix is given by Equation 
(73) and the stiffness matrix is given by 


m = 


0 i 0 

0 ' x m 


(78) 


and the [$f NN ] submatrix is defined in Equation (40). For the free-free 
mode synthesis the mass matrix is given by Equation (47) or (75) and the 
stiffness matrix by Equation (46). The damping matrix is a diagonal 
matrix of the form given by Equation (63). Equations (76) and (77) may be 
combined in the following uncoupled form: 


m { P }+ [O] f P }+ [X] { P } = [y(t)} 


(79) 


where 


m = 


TOo. 

(80a) 

[C] = 

>] B 

[<M 0- 

(80b) 

m = 


m 0 _ 

(80c) 

Cp] = \ 

<hh\ 
ltp3 0 / 

} 

(80d) 



(80e) 



Equation (79) now represents two matrix equations of motion for the 
unconnected Booster and Orbiter. The process of connecting the compo- 
nents by requiring that the displacements at their points of mutual attach- 
ment be equal is now discussed. 

2. Coupling Matrices . At the pin -joint attachment point we have 
the following displacement compatibility requirements 

x l,l ~ x i ,2 = 0 and 2 i,i ~ z i,b = 0 (81a) 

while at the slip-joint we have the displacement compatibility requirement 
that 


2 2,1 2 2,2 _ 0 


(81b) 


Using the pin- slip component modes and the component rigid-body modes, 
the x x>1 displacement may be written as follows: 


where 


-x,l 


= ClaPa 


+ <V,lrPr + 


m 

£ 

j=i 



(82a) 


0 


x 

1,1a 


normalized displacement of x 1jX in rigid-body 
axial mode 


0* >lr = normalized displacement of x lfl in rigid-body 

rotational mode 

0^^ j = normalized displacement of x t>l in the j-th 
elastic mode 

p°, p°, p° = generalized coordinates on Orbiter corresponding 
to the axial, rotational, and j-th elastic modes, 
respectively. 


Similarly, the remaining displacements may be written as 
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m 

*1,1 = + *LA + E ^IjP? < 82b > 

j = l 

n 

X l ,2 = ^l, 2 aPa + ^l*2 r Pr + E ^*l> 2 jpj (82 c) 

j = l 

*1- = *uajPl + + E ^4 «p! < 82d > 

j=l 

m 

z 2,l = *1,1 A + *3*rPJ + E 02,ljP? ( g2e ) 

j = l 

n 

Z 2*2 = ^2,2£p£ + ^2*2rPr + E ^2>2jPj (82f) 

j = l 


where 


= 

normalized displacement of z la in rigid-body 
lateral mode 

B 0 

Pi 'Pi = 

generalized coordinates on Booster and Orbiter 
corresponding to the rigid-body lateral modes 

B 0 

Pj’Pj = 

generalized coordinate on Booster corresponding 
to the j-th elastic mode 

etc. 


Note, however, that for the pin- slip case the elastic modal displacements 
at the boundary nodes were zero in the appropriate directions [Equations 
(31a, b) and (32a, b)]. Therefore, all the modal displacements in the 


summation in Equations (82a) through (82f) are zero since they correspond 
to these constrained displacements. From the compatibility requirements 
given by Equations (81a, b) we obtain three equations of displacement com- 
patibility in m + n + 6 generalized p coordinates where m = numberof 
Orbiter elastic modes and n = number of Booster elastic modes used in 
the analysis. These equations may be written in matrix form as 
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where 


[i: ' 


• 0 i 0* ! 0 ] 



0 




0 

“01 X ,2 r" 

[*!] = 

0 

01 l ,sJL 

"01 Z ,2r 


0 

02*2 

~0J,3 r_ 


<2> x 

W 1 ,1a 

0 

01 X ,lr _ 

[*!] = 

0 

$1 ,1 i 

0/ ,1 r 


0 

02,1 i 

02,1 f _ 


(83) 


(84a) 


(84b) 


and 


Pa’ 


[ p “i 

i 


Pa' 


0 

Pi 

1 

P 1 

> » C Pb n } H 

i 

1 

1 

l 

> > W =< 

A 

>. tp"3 •< 

1 

: 

i 

/r, 


p B 


Pr, 


i 

p ° 

m 


Due to the fact that there exist three equations of compatibility the com- 
ponents in the {p} vector are not independent. Recall that if there are 
& variables and there exists k equations of constraint between them, 
there will be i = 4-k independent variables or degrees-of-freedom. 

In our case we have m + n + 6 variables and three equations of constraint 
so that the number of degrees-of-freedom is m + n+3. It is entirely 
arbitrary as to which variables are selected as the dependent variables 
but in this study we choose the Booster generalized coordinates p , p. , 

B “Jo 

p r . Equation (83) may be written as 



from which 
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[0*] Cp*} + [A] {q} = 0 (86) 

where {q} is now a column vector of independent generalized coordi- 
nates. The dependent coordinates may now be obtained from Equation 
(86) as 


3 = "l>s I" 1 [A] Cq) 


(87) 


From this we derive an equation which relates the complete { p } vector 
to the independent £ q } sub-vector 


(p} = 



-Ck r 1 [a] 
to . 


Cq3 


( 88 ) 


The rectangular matrix in this equation is of order (m + n +6) X (m + n + 3) 
and [I] is an identity matrix of order m + n+3. Equation (88) may be 
written briefly as 


Cp} = [0] p Cq} (89) 

where [9] p is a coupling matrix analogous to the one defined by Hurty^ 
with the subscript p indicating that the pin-slip component modes were 
used for this analysis. 

For the synthesis of the free-free - component modes recall that 
all the displacements are referenced with respect to the interior nodes. 
Therefore, due to this and the fact that boundary nodes in the free-free 
component configurations have non-zero displacements we must first 
determine the boundary displacements in terms of the free-free rigid- 
body and elastic modes. Note that for this case the elastic components 
are not all zero as was the case for the pin- slip configurations. The 
displacement compatibility Equations (81a, b) still apply here but we must 
make the following additional requirements in order to transform the dis- 
placements at the interior nodes to the boundary nodes 

X x,x = x 3 + 6 i e 2 (90a) 

X 1 ,2 = x 8 - 6 1 0 8 (90b) 

Z X,1 = Z 2» Z X »2 = Z 8» Z 2,x = Z 6» Z 2,2 = Z XX (90c) 

The displacements x 2 , 9 S , x 8 , and 0 8 maybe written as 
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*3 = *;. a ° + 


0 3 = 0 3 , r P° + 


x 8 ~ $8* aPa + 


M 6 B 

0 8 = 0 8 ,rPr + 


m 

E 

*£jpj # 

j=i 


m 

r— * 

0 

z 


j=l 


n 


£ 


(1 

t— • 


n 

f— » 

6 a 

Z 

0 8,jPj 

j=i 



(91a) 

(91b) 

(91c) 

(91d) 


where 0 2(a tbe normalized displacement of x 2 due to rigid-body axial 
modes, etc. Similarly, the displacements given by Equation (90c) may be 
written as a combination of rigid-body and elastic modal displacements. 
Using Equations (91a, b, c, d) and Equations (90a, b, c) and substituting into 
Equations (81a, b) we arrive at a coupling matrix for the free-free compo- 
nent analysis which again relates the complete {p} vector to the indepen- 
dent £ q } sub-vector by 


£p3 = [P] f Cq)' (92) 

with the subscript f indicating that the free-free component modes are 
used in this analysis. Having obtained the coupling matrices which enable 
us to connect the Booster and Orbiter components to form the System we 
now proceed to the System eigenvalue problem. 

3. Eigenvalue Problem . As stated in the previous section, the 
coupling matrix [3] relates the complete £ p } vector to the independent 
sub-vector {q} which enables us to obtain the coupled System equations. 
Substituting either Equation (89) or Equation (92) (depending on whether 
we are synthesizing with the pin- slip component modes or the free-free 
component modes) into Equation (79), and then premultiplying by the 
transpose of the appropriate coupling matrix, we obtain the System 
equations of motion 


where 


t*l], t q } + [O], tq}+ [*], tq J = f Q(t) } s 


(93) 
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m 3 = m T m [3] 

(94a) 

[G], = [3] T [C] [3] 

(94b) 

r*]. = m T m m 

(94c) 

[Q(t)] s = [3 ] t {S(t)3 

(94d) 


Equation (93) is an (m + n+3) X (m+n+3) matrix equation where again 
m = number of Orbiter component modes and n = number of Booster compo- 
nent modes used in the analysis. To determine the frequencies and mode 
shapes for the System we solve [using Equation (93)] the undamped eigen- 
value problem 

i« 3 [I] - [D] =0 (95) 

where 

[d] = m ; 1 [*]. (96) 

A modified Jacobi method was used to obtain the eigenvalues and 
eigenvectors of Equation (95). The eigenvectors are normalized and may 
be written as an m + n+3 square matrix in the following form 

[Y] = [£Y x 3, {Y 2 3,....{Y m+n+3 3] (97) 

These eigenvectors may be used to relate the generalized coordinate 
vector { q } to a new vector { iq} which expresses displacements in a 
system of normal coordinates: 

Cq] = [Y] [ti] (98) 

At this point, however, we do not solve for the forced responses in the 
coordinate system. It is more useful to make a transformation back 
to our original coordinate system {u3. The reason for this transforma- 
tion at this point is twofold. First, it permits us to compare directly the 
analytically predicted mode shapes and forced responses with the re- 
sponses measured in the laboratory. Second, and most importantly, it 
permits us to apply a newly-developed method for synthesizing the compo- 
nent damping which appears to be extremely promising in terms of the 
accuracy obtained. This method of damping synthesis is discussed in 
detail in Section VI of this report. At this point we will only say that the 
method permits us to obtain System modal damping which results in a 
diagonal System damping matrix. This causes the equations of motion 
to be uncoupled and permits rapid solution of the forced responses. 
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For our 39 coordinate model (21 for the Booster and 18 for the 
Orbiter) we may obtain the System normalized mode shapes in the {u} 
coordinate system from the following equation 



\K ! o-i 



[0],„ [V] 


p ft 


(99) 


where again the subscript p denotes the pin- slip case and f denotes 
the free-free case. The System normalized modes are obtained by 
dividing each column of the resulting matrix product given in the above 
equation by an appropriate constant. The submatrices [0 B ] and [0 O ] 
are the component Booster and Orbiter modes used in the analysis, 
m is the coupling matrix, and [Y] is the eigenvector matrix given by 
Equation (97). The sizes of the matrices at this point are as follows: 

39 X (n + m + 3) 

39 X (n + m+ 6) 

(n + m+6) X (n + m+3) 

(n + m+3) X (n + m+3) 

Note that the System mode shapes [0] 3 are given in the {u} coordinate 
system which includes displacements at both interior and interface point 
positions. For convenience, we finally transform these mode shapes to 
the {u} or interior coordinates by means of Equation (21). That is, 
for the j-th column or mode shape we have 

U,} s = [a] {*,}, 

The displacements (x lf z t , 9j) at node i may now be expressed as a 
summation of system modes and system generalized coordinates, i. e. , 



n+m 


X 1 = <0! + K iVl + 0u T Pr + £ ^kPk 

k= 1 


n+m 


z ! = 0l\aPa + 0\, m + ^i.rPr + £ ^l.kPk 

k=l 


n+m 


0 j 9 3 j9 S ^ 9 S V rk ® S 
9 1 = ^l,aPa + *i,m + 0 i,rPr + £ ^i,kPk 

k=l 


(100a) 


(100b) 


(100c) 
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or letting 


(u^s = 


( 101 ) 


we may write Equations (100a, b, c) in matrix form as 

0^] s 



( 102 ) 


For all the nodes on both the Booster and the Orbiter we may 
combine and arrange the displacements in the following manner 


{u} s = 


tu} 0i 


(103) 


where 



( _ ' 





u 7 



U 1 

{u} 6 = < 

U 8 

1 

1 

^ and 

II 

O 

V 


i 



U 6 


(104) 


which are synonomous with Equation (8). Likewise, Equation (102) may 
be written to include all the nodes as follows: 


{uL = 


‘k I 
1 

0 N " 
w B 

jl > 

Z N 

0 oJ 


= Ms tph 


(105) 


where the first 21 rows of the modal matrix correspond to Booster modal 
diaplacements and the next 18 rows correspond to Orbiter modal displace- 
ments with each mode in the matrix arranged column-wise. The general- 
ized System mass and stiffness matrices may now be written as 


m s = 

~K 

i 

1 

j_ 

K 

T 

”m b 

1 

1 

. L 

0 ~ 



i 

i 

i -o-i 
i — , 

?o 

1 

i 

1 

K 


0 

i 

i 

i 

M 0 _ 



i 

it 

CO 

*2 


1 

1 

l 


T 

K. 

t 

i 

0 


~~K 


K 

T 

\ 

i 

0 N 
^ °_ 


0 

i 

i 

« 

K 0 . 


jl 

1 V oj 


(106a) 


(106b) 
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where the tilde (~') indicates that these quantities are with respect to 
the System normal modes. 

4. Solution of Coupled Equations . Equation (93) represents the 
System equations of motion in matrix form and at this point is highly 
coupled. Solution of this system of coupled, damped equations would 
appear to be a straightforward possible method of damping prediction 
for the System. It was decided to pursue this possibility by solving for 
System modal damping. A check on the damping matrix was first per- 
formed by solving the damped eigenvalue problem. A solution to the 
undamped eigenvalue problem showed very good agreement with the 
natural frequencies and mode shapes obtained on the model in the labora- 
tory. A method for solving equations of motion for a damped linear system 
was developed by Foss(^). The procedure is fairly lengthy and only the 
salient features of the method will be discussed here. 

For a damped system having n equations of motion, the eigen- 
value and eigenvectors are complex for a system that is less than criti- 
cally damped. This means that 2n equations are required since two 
items of information (magnitude and phase angle) are required to deter- 
mine each mode. Following Foss' method we write the following matrix 
identity 


m, {q} - [»!], tqJ - (0} (107) 

where [^] s is the mass matrix shown in Equation (93). This equation 
is combined with Equation (93) which is repeated here 

M, tq} + [<J] s Cq} + [>Os fq) = Cat*)}, ( 93) 


The two equations are then combined as follows: 

[S] {y} + [T] {y } = { Y(t) } 


where 


[S] 

[T] 


£y3 


" o 1 m. 

i_ _ _ 

», ! » s 

L. 

X ! 0 ' 

_° S *r s 



(108) 


(109a) 


(109b) 


(109c) 



38 


Cy (t)} = 



(109d) 


The solution to the homogeneous part of Equation (108) is obtained by 
setting the right-hand side equal to zero and may be expressed in the 
exponential form, e pt . Therefore, we have the following relationship: 


{y 3 = p {y 3 


which when substituted into Equation (108) gives 

p[S] (y 3 = — [T] Cy} 


( 110 ) 


dll) 


Equation (111) is an eigenvalue problem that yields a set of 2n eigen- 
values and a set of 2n eigenvectors. The imaginary parts of the eigen- 
values give the frequencies of the free vibrations and the real parts give 
a measure of the rate of exponential decay of the amplitudes of free 
vibration. The eigenvalues may be written as follows: 

P, = -C, »! 1 (112) 

where id* is the damped natural frequency for the j-th mode and Cj is 
the damping factor for the j-th mode which can be determined by dividing 
the real part of p^ by ID*'. It should be noted that Gj, as determined from 
the complex eigenvalue procedure, is identical to that which should be expected 
to result if the coupled forced response were solved directly, and Gj were 
determined from a g-bandwidth procedure applied at each resonant frequency. 


Comparison of the damping factors obtained by solving the above 
eigenvalue problem with those obtained on the model showed the results 
to be extremely poor with the majority of errors being in the neighborhood 
of 100-325 percent, as will be shown in Section VII. B. As a result of the 
poor agreement, solution of the nonhomogeneous equations would have pro- 
duced equally poor results especially at the System natural frequencies, 
since here the responses are damping controlled. Anticipating this fact 
at the outset of this program, we studied various other ways by which the 
System damping could be synthesized from the component damping with 
favorable results. A new method for damping synthesis was developed 
which produced extremely good results, at least for our present fairly 
simple model. Details of this method will be given in Section VI. How- 
ever, we first conclude the theoretical development by indicating how 
predicted modal damping is to be used for predicting forced response. 
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5. System Forced Response Using Predicted Modal Damping. 

Using the modal mass and stiffness matrices^given by Equations (106a, b), 
we now introduce a modal damping matrix [C-] 3 which is obtained from 
the procedures outlined in Section VI of this report. This damping matrix 

is diagonal with the j-th element given by C' i = 27l\ j cjUj Cj where is 
the j-th natural frequency obtained from the eigenvalue problem given 
by Equation (95) and is the predicted damping factor for the j-th 
System mode. The uncoupled equations of motion are now written in 
matrix form as 

m, CpU [«], ({>},+ [?], (p)„= t Q(t) 3 , (113) 

where {Q(t)} s is a vector of generalized applied forces. Procedures 
for obtaining the System forced responses are the same as for the compo- 
nent forced responses which were given in Section III. C. It should be noted 
that the damping factors used in the System forced response analyses are 
the predicted ones and not those obtained experimentally, although a com- 
parison of the two is given in a later section. Once the generalized {p} 
coordinates are known, the displacements may be obtained by substituting 
into Equation (105). 



IV. EXPERIMENTAL PROCEDURES 


A. Background Discussion 

In order to utilize damping results from component tests for 
prediction of System behavior, it was recognized that some form of 
correlation of component results would first have to be established. 
Several methods of doing so were investigated. The method finally 
found to work best is based on a modification of a damping energy 
correlation which was presented in Reference 1. However, we now 
will show that the present modification is necessary in order for the 
correlation to be valid. 


A straightforward dimensional analysis of the component damp- 
ing results can be performed by considering the following variables: 

D o = Damping energy dissipated per cycle 

M = Modal mass 

uu = Natural frequency of a given mode 

x q = Peak modal amplitude for steady- state sinusoidal 
vibration 

X = Some characteristic length 


A dimensional analysis results in the nondimensional equation 


D„ 


M<JJ 3 x 


J—S = ^ x {x 0 /i) or 


D„ 


= (* 0 / 4 ) 


where T 0 = ijMU) 3 x 0 2 is the peak steady state kinetic energy at ampli- 
tude x 0 . This equation can be written in dimensional form as 


D e = * a (T o, *«,/*) (1 

This result shows that damping energy is a function of two independent 
variables, not just kinetic energy alone, as was originally assumed in 
Reference 1. 


In view of Equation (114), single-curve correlation will result 
only if x 0 /£ is held constant, and a family of curves is obtained, one 
for each amplitude. Sample results for one Booster configuration are 



shown in Figure 5, plotted in the form of Equation (114). Strictly- 
speaking, the most general procedure for obtaining. such data is to 
measure input energy at each of the component resonances, while 
holding the amplitude constant. Kinetic energy is computed from 
measured frequencies and theoretical modal mass. However, for 
a system with many resonances, it is extremely difficult to obtain 
accurate measurements at low frequencies, for those amplitudes 
which are appropriate at higher modes. On the other hand, if the 
damping forces are found to vary essentially linearly with ampli- 
tude x 0 , and damping energy linearly with x 0 s , then measurements 
can be made at any convenient amplitude, and adjusted to the values 
corresponding to a given value of x o /i. Further, if this condition 
does exist, then it is also possible to measure a modal damping ratio 
C by means of conventional free decay or half -bandwidth techniques, 
and the damping energy can be calculated from 

D c = 4ttT 0 £ 

The latter procedure was utilized in reducing the majority of the data 
in the present program. However, it should be noted emphatically, 
that data such as that shown in Figure 5 can be generated for any form 
of damping process (including one nonlinear with amplitude), so long 
as damping energy is measured directly, and amplitudes are either 
held constant, or values are extrapolated by some suitable procedure. 

There is another extremely important point which must be 
recognized with regard to the form of plot displayed in Figure 5. That 
is, the points represent measurements made only at resonant conditions 
for a given component. When a curve is drawn through these points, 
we automatically imply that damping energy can be obtained at inter- 
mediate, off-resonant kinetic energies as well from these curves. In 
other words, since we have a system with light damping and reasonably 
separated modes. Figure 5 is not only a plot of modal damping energy 
versus peak modal kinetic energy, but it can be hypothesized to be a 
plot of total damping energy versus total peak kinetic energy as well, 
since at resonance, only one mode dominates. The latter statement 
must hold true regardless of the vibrational form of the component. 

It would be appropriate to consider measurements at off- resonant con- 
ditions to obtain further data to support this assertion, since it will be 
a fundamental hypothesis to be used in the later development of our 
method for prediction of System modal damping. However, measure- 
ments at such conditions are usually highly inaccurate, and other further 
evidence which confirms this hypothesis will be given near the end of 
the report. 








B. 


General Procedures 


It is recognized that the primary emphasis of this program is to 
develop a method for prediction of System damping. However, in order 
to do so, use of other dynamic modal properties of the components and 
the System must be included. It is logical that natural frequencies and 
mode shapes must be considered as part of the comparison of results 
from analytical and experimental models. Likewise, comparison of 
results for steady-state forced response are also desirable in addition 
to values for modal damping. Therefore, we will include data to show 
the effects of a variety of parameters on component behavior, and in 
many cases will also include predicted values for the components for 
comparison. These results will establish a confidence in the accuracy 
of the component mathematical models. We subsequently will then dis- 
cuss results for the System. 

The approximate order of conducting the experiments has been 
indicated in Table 1(a). Once the apparatus for a given configuration was 
set up, natural frequencies and mode shapes were measured first. Exci- 
tation was provided by a very light electrodynamic exciter which was 
located at Mass 10 on the Booster and at Mass 4 on the Orbiter. For 
System configurations part of the data were obtained by driving on each 
component. Mode shapes were measured in the lateral directions only 
by means of multiple piezoelectric accelerometers which were mounted 
to various masses. Most of the accelerometer cables were removed for 
subsequent damping measurements, since it was found that motion in 
these cables could cause damping of the same order as the very light 
structural damping exhibited by the components when tested in air. 

Damping measurements were performed by the 1 /2-bandwidth 
technique at first. Only slight nonlinearity with amplitude was usually 
encountered. However, those configurations involving the lightest 
damping could best be tested using the free decay method. Most of the 
data were subsequently obtained using this method, since a permanent 
record resulted, and the results were commensurate with those obtained 
from 1 / 2 -bandwidth measurements. 

During this entire program, careful measures were taken to en- 
sure the reliability of the results. Oil temperatures and viscosities were 
monitored several times daily. Independent calibration of the force 
exciter and accelerometers were repeated on several occasions. Correc- 
tions to damping values to allow for viscosity changes were based on 
day-to-day plots of the data. Tests on various modes were repeated at 
later times when possible in order to provide cross-checks in the data. 
Sample data were often taken at various amplitudes to spot any significant 
nonlinearity. After all these measures, it is estimated that the worst 
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error for modal damping measurements is ±15%, although most of the 
data is much less. This is not bad, however, as damping measurements 
on typical structures are usually much worse. It is also surmized that 
only the very first few test configurations experienced much of the worst 
error, as improvements in procedures were incorporated as time pro- 
gressed and behavior of the apparatus was better established. 

Tests for properties in rigid-body modes were necessary to pro- 
vide data for System damping prediction. This was a consequence solely 
of our particular design of providing damping; that is, damping relative 
to the fixed ground. Axial and lateral tests were conducted by measuring 
decay rates at the low frequencies (about 1 Hz) of the models on the 
simulated free-free suspension system. Rotational tests were conducted 
at about the same frequency by supporting the models on a knife-edge at 
their CG, and stabilizing the model with a small lateral spring. In all 
cases, the same oil levels were maintained as for the tests for elastic 
modes. 

C. Link Damper Calibration 

The link damping mechanism described in Figure 3 evolved from 
a series of other mechanisms which proved to be far more unreliable. 

It would have been desirable to have calibrated this device with a com- 
pletely independent apparatus. However, in order to save time, it was 
decided to calibrate it directly in place on the combined System. It must 
be emphasized that this in no way deters from the validity of the results 
that followed. 

For calibration, the link damper was installed on the combined 
System, while it was supported in air (i. e. , in the DO condition with no 
oil in the reservoirs). Relative velocity across the link was monitored 
at a steady- state forced response for each mode, first in the LI condition 
(with no feedback in the damper). Measurement of the input energy thus 
allows establishing a small damping energy reference level associated 
with the internal structural stresses and the motion of the flexure pins. 

This level was usually less than 0. 1% critical in any mode. Then, the 
energy level with feedback on was also measured at the same input 
amplitude. The difference of the two energy levels was that added by 
the link damper, and was reduced to an equivalent viscous damping 
coefficient, since the measured relative velocity across the link was 
available in all cases. This procedure resulted in the damping coeffi- 
cients listed in Table III. A further check during these procedures showed 
that no significant nonlinearity or change in mode shapes occurred for the 
various levels of damping employed. This was somewhat surprising, since 
at the largest damping level (L4), the joint provided more than 30 times 
the damping energy absorbed in the rest of the structure when it was 
supported in air. 
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TABLE III. LINK DAMPING COEFFICIENTS FOR L2, L3, L4 


Mode No. 

C U) 

C (a) 

c (3) 

1 

93. 05 

185. 83 

354.63 

2 

108. 80 

212. 78 

523.43 

3 

122. 98 

235.47 

434. 07 

4 

77. 59 

161. 71 

374.49 

5 

123. 13 

270. 94 

509. 25 

6 

113. 20 

178. 73 

431. 23 

7 




8 



___ 

9 

91. 92 

270. 94 

187. 24 

10 

_ _ — 

_ _ _ 

_ 


(1) 

L2 

Link Damping 

Units for C are newton-sec/meter 

(2) 

L3 

Link Damping 



(3) L4 


Link Damping 
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V. COMPONENT RESULTS 


A. Pin -Slip Configuration 

1. Natural Modes . The experimental and predicted normalized 
mode shapes for the Booster in the Ml and M2 mass configurations are 
shown in Figures 6a and 6b, respectively. Only the lateral direction is 
presented although longitudinal and rotational motions were also present 
through coupling. Six mode shapes were determined experimentally for 
each mass configuration. For the Orbiter in the Ml and M2 mass con- 
figurations, the normalized mode shapes are shown in Figures 7a and 7b, 
respectively. All the experimental mode shapes are shown in the D2 
damping condition. Due to the slight difference in densities of the various 
oils used to produce the damping, it was surmised that there might be a 
change in mode shapes between damping conditions. Comparison of the 
mode shapes among the various damping conditions showed, however, that 
any changes were insignificant. The experimental and theoretical natural 
frequencies and modal masses are presented in Table IV(a) for the Booster 
and in Table IV(b) for the Orbiter. Results are presented for various 
damping conditions to show the effect of the various oil densities on the 
natural frequencies of the component. The theoretical frequencies given 
were obtained for the D2 oil density condition. 

It was discovered during the course of experimentation that there 
was a small amount of translational compliance at the pin-joints. This 
seemed to affect particularly both the Booster and Orbiter fourth modes, 
with the difference in theoretical and experimental frequencies being about 
7 to 8 percent. The percent difference is based on experimental fre- 
quencies. This difference was not considered exceptionally bad, and it 
was further surmised that this discrepancy would not influence damping 
results. Therefore, no modifications were made to the pin-joints. 

2. Damping . Experimental damping energy per cycle versus 
kinetic energy curves are shown in Figures 8a and 8b for the Booster 
and Orbiter substructures, respectively. The curves were obtained by 
fairing through the experimental data points given by the symbols, as 
has already been discussed in Section IV. A. A least squares fit would 
have been appropriate but was not considered here. Results are shown 
for the DO, D2, and D3 damping conditions in the Ml mass configuration 
and for the D2 damping condition in the M2 mass configuration. It 
should be noted that both the Ml and M2 mass configuration energy curves 
appear to be the same (at least within the scatter of the data) for the D2 
condition. This result was particularly encouraging since it seems to 
indicate that for different mass -configured substructures with similar 
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Figure 7. Orbiter Mode Shapes for Pin- Slip Configurations 

(a) Ml Mass Condition 




TABLE IV (a). COMPONENT NATURAL FREQUENCIES AND 
MODAL MASS FOR PIN- SLIP SUPPORT 

( Booster ) 


Configuration 

Mode Number 

B-P-L1-D1-M1 

1 

2 

3 

4 

5 

6 

■j* 

Theoretical frequency (Hz) 

6. 18 

13. 12 

33. 09 

50. 68 

59.96 

94. 14 

Experimental frequency (Hz) 

6. 31 

13. 12 

33.43 

47. 00 

58. 00 

91.40 

Percent difference 

-2. 06 

0. 00 

-1. 02 

+ 7. 83 

+ 3. 38 

+ 3. 00 


B-P-L1-D2-M1 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 

6. 18 

13. 12 

33. 09 

50. 68 

59.96 

94. 14 

Experimental frequency (Hz) 

6. 26 

13. 09 

33. 20 

46. 90 

57. 80 

90. 90 

Percent difference 

-1.28 

+0.23 

-0. 33 

+8. 06 

+ 3.74 

+ 3. 56 


B-P-Ll -D3-M1 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 

6. 18 

13. 12 

33. 09 

50. 68 

59.96 

94. 14 

Experimental frequency (Hz) 

6. 20 

13. 06 

33. 15 

46. 80 

57.90 

91. 00 

Percent difference 

O 

CO 

• 

o 

1 

+ 0. 46 

-0. 18 

+8. 29 

+ 3.56 

+ 3. 45 


B-P-L1-D2-M2 

1 

2 

3 

4 

5 

6 

•j* 

Theoretical frequency (Hz) 

6. 57 

13. 28 

35.48 

52. 37 

69. 09 

150. 06 

Experimental frequency (Hz) 

6. 64 

13. 75 

35. 62 

48. 20 

66. 20 

122. 90 

Percent difference 

-1. 05 

-3. 42 

-0. 39 

+ 8. 65 

+4. 36 

+22. 10 


Normalized Modal Mass 

1 

2 

3 

4 

5 

6 

Ml mass configuration (Kgs) 

2. 50 

3. 55 

4. 85 

9. 60 

4. 32 

3. 52 

M2 mass configuration (Kgs) 

2. 21 

3.47 

3. 15 

12. 96 

2. 73 

2. 78 


v 

Undamped eigenvalue results 








TABLE IV (b). COMPONENT NATURAL FREQUENCIES AND 
MODAL MASS FOR PIN- SLIP SUPPORT 

( Orbiter ) 


Configuration 

Mode Number 

O-P-Ll-Dl-Ml 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 

10. 52 

18. 58 

37. 22 

57. 34 

93.86 

_ _ _ 

Experimental frequency (Hz) 

10. 62 

18. 40 

36. 73 

53. 30 

92. 10 

_ _ _ 

Percent difference 

-0. 94 

+ 0. 98 

+ 1. 33 

+7. 58 

+ 1. 91 




0-P-L1-D2-M1 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 

10. 52 

18. 58 

37. 22 

57. 34 

93. 86 

_ _ _ 

Experimental frequency (Hz) 

10. 67 

18. 35 

36. 67 

53. 10 

92. 00 

_ _ _ 

Percent difference 

o 

• 

1 

+ 1. 25 

+ 1. 50 

+7. 98 

+2. 02 




O-P-Ll - D3-M1 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 

10. 52 

18. 58 

37. 22 

57. 34 

93.86 

_ _ _ 

Experimental frequency (Hz) 

10. 58 

18. 30 

36. 62 

53. 10 

92. 20 

_ _ _ 

Percent difference 

-0.57 

+ 1. 53 

+ 1. 64 

+ 7. 98 

+ 1. 80 

— 


0-P-L1-D2-M2 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 

10. 71 

18. 61 

39. 28 

57. 07 

114. 05 

_ _ _ 

Experimental frequency (Hz) 

10. 73 

18. 40 

38. 55 

53. 60 

106. 80 

_ _ _ 

Percent difference 

-0. 19 

+ 1. 14 

+ 1.89 

+6.47 

+6. 69 

— 


Normalized Modal Mass 

1 

2 

3 

4 

5 

6 

Ml mass configuration (Kgs) 

5. 57 

3. 17 

5. 99 

8. 90 

3. 38 



M2 mass configuration (Kgs) 

5. 29 

3. 06 

3. 96 

6. 71 

2. 61 

— 


Undamped eigenvalue results 
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type damping mechanisms, the damping energy curves are congruent. 
Modal damping ratios in terms of critical damping are given 
in Table V for both the Booster and Orbiter for various mass and 
damping conditions. Values for the DO (in air) conditions are also given 
to indicate that most of the damping was provided by the discrete oil 
dampers. The DO values for the Orbiter tend to be somewhat high, 
and are attributed to extraneous damping of accelerometer cables which 
was eliminated in subsequent tests. 

3. Forced Response . To obtain forced responses, the Booster 
was excited with a constant force at Mass 10 in the lateral (z) direction 
and the Orbiter at Mass 4. Response measurements were also made at 
these locations in the lateral direction by means of an accelerometer 
attached to the component. Figures 9a and 9b show the experimental 
and predicted forced responses for the Booster and Orbiter, respectively, 
in the form of response in DB versus frequency. The responses are 
given in terms of velocity per unit force, with the response in the first 
mode taken as the reference. Therefore, 

V/F 

DB = 20 log 

where V is the velocity at frequency f and V x is the velocity in the 
first mode, while the force F is held constant. These forced responses 
serve as an indication of the accuracy of-the analytical prediction proce- 
dure. 

B. Free-Free Configuration 

1. Natural Modes . The experimental and predicted normalized 
mode shapes for the Booster and Orbiter in the Ml mass configuration 
are shown in Figures 10 and 11, respectively. Here again, as for the 
pin- slip configuration, these mode shapes are presented for the D2 damp- 
ing condition. The experimental and theoretical frequencies and modal 
masses are presented in Table VI. Included are the three frequencies 
associated with the suspension system natural frequencies and denoted by 
A = axial, L = lateral, and R = rotational rigid-body motions, respectively. 
As noted from Table VI, agreement between experiment and theory is very 
good for the six Booster and five Orbiter modes. The percent difference 
between experimental and predicted natural frequencies is again based on 
the experimentally measured frequencies. 

2. Damping. Experimental damping energy per cycle versus 
kinetic energy curves are shown in Figures 12a and 12b for the Booster 
and Orbiter substructures, respectively. The curves were obtained by 
fairing through the experimental data points given by the symbols, and 



TABLE V. COMPONENT DAMPING FOR PIN- SLIP SUPPORT 

( Booster and Orbiter ) 
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Configuration 

Mode Number 

B-P-Ll -DO-MI 

I 

2 

3 

4 

5 

6 

. 0008 

. 0010 

. 0011 

. 0012 

. 0014 

. 0016 

O-P-Ll-DO-Ml 

. 0024 

. 0016 

. 0032 

. 0024 

. 0020 




B-P-Ll -D1 - Ml 

. 0092 

. 0057 

. 0047 

. 0034 

. 0046 

. 0032 

O-P-Ll -D1 -Ml 

. 0079 

. 0085 

. 0074 

. 0038 

. 0045 

- 


B-P-L1-D2-M1 

. 0164 

. 0101 

. 0066 

. 0056 

. 0087 

. 0051 

Q-P-L1-D2-M1 

H9 

. 0117 

. 0096 

. 0053 

. 0058 




B-P-Ll -D3-M1 

. 0365 

. 0200 

. 0119 

. 0071 

. 0191 

. 0080 

Q-P-L1-D3-M1 

. 0244 

. 0223 

. 0142 

. 0070 

. 0087 




B-P-Ll -DO -M2 

. 0013 

. 0025 

. 0019 

. 0026 

. 0022 

. 0033 

O-P- LI - DO -M2 

. 0020 

. 0015 

. 0037 

. 0028 

. 0017 

— 


B-P-L1-D2-M2 

. 0158 

. 0112 

. 0074 

. 0063 

. 0101 

. 0116 

Q-P-L1-D2-M2 

.0121 

. 0128 

. 0086 

. 0047 

. 0075 

— 
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Figure 9a. Booster Component Response at Mass 10 
for B-P-Lil -D2-M1 Configuration 
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Figure 9b. Orbiter Component Response at Mass 4 
for 0-P-L1-D2-M1 Configuration 











MODE 3 MODE 4 


Figure 11. Orbiter Mode Shapes for Free-Free Configurations 

Ml Mass Condition 


TABLE VI. COMPONENT NATURAL FREQUENCIES AND 
MODAL MASS FOR FREE- FREE CONDITION 

( Booster and Orbiter ) 


Configuration 

Mode Number 

B-F-L0-D2-M1 

A 

L 

R 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 



_ _ _ 



8. 53 

24. 08 

45. 72 

75. 25 

106. 10 

203. 40 

Experimental frequency (Hz) 

0.60 

3. 15 

3. 16 

8. 66 

24. 41 

46. 20 

76. 20 

109. 20 

196. 10 

Percent difference 

— 





-1. 50 

-1. 35 

-1. 04 

-1. 25 

-2. 84 

+ 3. 72 


O-F-L0-D2-M1 

A 

L 

R 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 



_ .» _ 



11. 61 

33. 11 

64. 56 

98. 81 

203. 90 



Experimental frequency (Hz) 

0.77 

3. 57 

3. 75 

11. 72 

, 

33. 19 

64. 60 

100. 70 

200. 10 



Percent difference 







-0. 94 

-0. 24 

-0. 06 

00 

00 

+ 1. 90 

— 


B-F-L0-D3-M1 

A 

L 

R 

T1 

2 

3 

4 

5 

6 

❖ 

Theoretical frequency (Hz) 







8. 53 

24. 08 

45. 72 

75. 25 

106. 10 

203.40 

Experimental frequency (Hz) 

0.60 

3. 15 

3. 16 

8. 62 

24. 50 

46. 10 

76. 30 

109. 10 

196. 20 

Percent difference 

— 

— 

— 

-1. 04 

-1.71 

-0. 82 

-1. 38 

-2. 75 

+ 3.67 


O-F-L0-D3-M1 

A 

L 

R 

1 

2 

3 

4 

5 

6 

Theoretical frequency (Hz) 

. _ _ 





11. 61 

33. 11 

64. 56 

98. 81 

203. 90 



Experimental frequency (Hz) 

0.72 

3. 58 

3. 23 

11.68 

33. 20 

64. 60 

101. 00 

201. 50 



Percent difference 

— 


— 

-0. 60 

-0. 27 

-0. 06 

-2. 17 

+ 1. 19 

— 


Normalized Modal Mass 

A 

L 

Rt 

1 

2 

3 

4 

5 

6 

Booster 

Ml mass configuration (Kgs) 

12. 26 

12. 26 

4. 74 

6. 46 

6. 65 

6. 30 

7. 22 

4. 69 

9. 18 

Orbiter 

Ml mass configuration (Kgs) 

10. 44 

10. 44 

5. 72 

6.62 

7. 37 

4. 69 

5. 71 

10. 24 




A = Axial rigid body L = Lateral rigid body R = Rotational rigid body 

Undamped eigenvalue results 
^ If l 2 (l - distance from eg to bottom of component) 





















results are shown for the DO, D2, and D3 damping conditions in the Ml 
mass configuration. In additon to the damping values associated with 
the natural modes, damping values are plotted for the rigid-body motions 
associated with the axial, lateral, and rotational suspension system 
natural frequencies. Recall that due to the nature of the damping 
mechanism employed in the experiments, damping in rigid-body modes 
is also present at the non-resonant frequencies of the component. It 
was then assumed for the purpose of analysis that the damping associated 
with rigid-body motion corresponded to the damping obtained at the sus- 
pension system natural frequencies. Experimentally determined elastic 
mode damping ratios as well as rigid-body damping ratios are given 
in Table VII for the Booster and Orbiter substructures. 


3. Forced Response. To obtain forced responses, the Booster 
was excited with a constant force at Mass 10 in the lateral (z) direction 
and the Orbiter at Mass 4. Response measurements were also made at 
these locations in the lateral direction by means of an accelerometer 
attached to the component. Figures 13a and 13b show the experimental 
and predicted forced responses for the Booster and Orbiter, respectively, 
in the form of response in DB versus frequency. As for the pin- slip 
configuration, the responses are given in terms of velocity per unit force, 
with the response in the first elastic mode taken as the reference value. 
As with the pin- slip configurations, these responses serve only as an 
indication of the accuracy of the analytical model, and offer nothing new 
in dynamic analysis. New findings will be presented for the System 
results. 
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VI. SYNTHESIS OF SYSTEM DAMPING 


The comments in Section IV. A provide the basis for our method 
of damping synthesis. In particular, we hypothesize that the curves in 
Figures 8a, b and 12a, b provide a means of selecting an appropriate value 
of damping energy for each component, once we establish its peak kinetic 
energy of vibration at a given frequency and amplitude, regardless of the 
spatial shape it may sustain. In view of this hypothesis, we are now in 
position to predict the damping values for the combined System. We 
first perform an ordinary modal analysis of the undamped System. This 
provides the natural frequencies, mode shapes, and kinetic energies for 
a given modal amplitude x 0 . Let T oS , T oB , and T o0 be the kinetic 
energies for the System, Booster, and Orbiter, respectively, in a given 
System mode. We note that either the Booster or Orbiter will have an 
amplitude x o while the other component will have some smaller maxi- 
mum amplitude in that mode. Since the energy curves (Figures 8a, b and 
12a, b) are plotted with constant amplitude x Q , the peak displacement of 
the component having the smaller displacement is scaled so that its peak 
displacement is now x o . (This is equivalent to selecting the damping 
energy from the curve of lower constant amplitude, were such a curve 
available. ) The scale factor is obtained by dividing x o by x' where 
x' <x q , and is the maximum amplitude in a given mode on the component 
having the smaller amplitude. The scale factor may then be written 

S f = x 0 /x' > 1 

For the component having amplitude x 0 we enter the appropriate energy 
curve and read the damping energy, D 0 , corresponding to the kinetic energy 
T 0 = M(ju 2 x 0 2 . For the component having an amplitude x' we first multiply 
the kinetic energy by the scale factor squared so that T 0 ' = i?Mou s x / 2 S f 2 . 

We then enter the appropriate energy curve and read a damping energy, D 
This damping energy, however, is based on an amplitude x . To obtain 
the actual damping energy we simply divide D e by the scale factor squared 
to obtain actual damping energy, D c . Once the D 0 values are determined 
for each component the System damping can then be predicted from 


C. 


^cB ^eO 

T 0$ 


c S 


4tt T 


o s 


(115) 


The damping ratios obtained using this method are presented in Section 
VII. B for several system configurations. It should also be pointed out that 
the above scaling procedure is necessary, since Figures 8 and 12 indicate 
that damping energy is a nonlinear function of peak kinetic energy. 
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Allowance for joint damping is accomplished by adding the damping 
energy of the joint which is calculated from the calibrated properties of the 
joint damper (see Section IV. C) and the relative velocity across the joint 
as given by the theoretical mode shape. The damping ratio for both the 
link and structure is then obtained by 


c 


8 9 


i 


D es + D =i 


C 6 + C i 


4tt T 


( 116 ) 



VII. SYSTEM RESULTS 


A. Natural Modes 


System natural frequencies and mode shapes were obtained from 
synthesis of both the component pin- slip and component free-free modes. 
For the Ml and M2 mass configurations the first six pin- slip Booster 
and first five pin- slip Orbiter modes were used in the synthesis. The 
first five Booster free-free modes and the first four Orbiter free-free 
modes were used in the synthesis for the Ml mass configuration only. 
Inspection of Ml mass configuration System modes and frequencies 
showed that the pin- slip component mode synthesis produced much better 
results at higher frequencies than did the free -free component mode 
synthesis. Therefore, only the pin- slip component modes were used in 
the M2 mass configuration System synthesis. 

The System mode shapes obtained from pin- slip component modes 
and free-free component modes are shown in Figure 14a for the Ml mass 
configuration. It can be seen that the agreement between the experimental 
and predicted values is very good to excellent for the pin- slip case up to 
the tenth mode. It should be noted that the seventh and eighth System mode 
could not be determined experimentally, because of the location of the 
exciter. However, since the pin-slip component modes produced good re- 
sults up to the tenth mode it was surmised that modes 7 and 8 would con- 
form to the pin- slip analysis and are so indicated. Figure 14b shows the 
mode shapes obtained for the M2 mass configuration, and they also agree 
quite well with the experimental results. 

Table VIII(a) shows the System natural frequencies and modal 
masses obtained from synthesis of the pin- slip component modes for the 
Ml and M2 mass configurations. It should be noted that three eigen- 
value results are presented for the Ml mass configuration. These re- 
sults correspond to 1) an analysis assuming zero damping, 2) an 
analysis using the D2 component damping results, and 3) an analysis 
using the D3 component damping results. These results were obtained 
from a complex eigenvalue analysis in which we tried to determine if 
System damping could be obtained by synthesizing the component damping 
matrices as outlined in Section III. D. 4. The differences in predicted 
frequencies for the undamped and the damped cases are, of course, very 
small and are presented here only for comparison. To assess the con- 
vergence properties inherent in assuming a finite number of component 
modes for synthesis, two other cases are presented in Table VIII(a), in 
which we synthesized three Booster and three Orbiter modes and then 
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TABLE VIII(a) 

SYSTEM NATURAL FREQUENCIES AND MODAL MASS 
FROM SYNTHESIS OF PIN- SLIP COMPONENT MODES 


Ml Mass Configuration 


Booster Modes 1 through 6 used. 
Orbiter Modes 1 through 5 used. 

Mode Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Theoretical frequency* (Hz) 

7.31 

J 1 . 02 

16. 25 

27. 38 

35. 74 

49. 63 

55. 87 

74. 92 

85. 67 

108. 81 

Theoretical frequency^ (Hz) 

7. 30 

11 . 01 

16. 25 

27. 37 

35. 73 

49. 63 

55.86 

74. 91 

85. 67 

108. 79 

Theoretical frequency"* (Hz) 

7.29 

11.01 

16. 24 

27. 33 

35.69 

49.62 

55. 83 

74.91 

85. 65 

108. 72 

Experimental frequency (Hz) 

7. 35 

11. 13 

16. 34 

26. 79 

34. 56 

49. 73 

56. 00 

76.00 

87.07 

108. 64 

4 

Percent difference 

-0. 54 

-0.99 

-0. 55 

+ 2 . 20 

+ 3. 41 

- 0 . 20 

-0.23 

-1. 42 

-1.61 

+ 0 . 16 


Ml Mass Configuration 


Booster Modes 1 through 3 used. 
Orbiter Modes 1 through 3 used. 

Mode Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Theoretical frequency* (Hz) 

7. 31 

11 . 02 

16. 37 

27. 78 

36.06 

55. 05 





Experimental frequency (Hz) 

7. 35 

11. 13 

16. 34 

26. 79 

34. 56 

49. 73 









A 

Percent difference 

-0. 54 

-0. 99 

+ 0 . 18 

+ 3.69 

+4. 34 

+ 10. 70 










Ml Mass Configuration 


Booster Modes 1 through 2 used. 
Orbiter Mode 1 used. 

Mode Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Theoretical frequency* (Hz) 

7. 65 

11 . 20 

16. 37 








Experimental frequency (Hz) 

7. 35 

11. 13 

16. 34 















4 

Percent difference 

+4. 08 

+0.63 

+ 0 . 18 





— 

— 

— 

— 

— 


M2 Mass Configuration 


Booster Modes 1 through 6 used. 
Orbiter Modes 1 through 5 used. 

Mode Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Theoretical frequency* (Hz) 

7.68 

11. 23 

16. 22 

26. 41 

38. 30 

50. 18 

56. 35 

88 . 53 

96. 89 

161. 65 

Experimental frequency (Hz) 

7. 75 

11. 33 

16. 48 

26. 19 

38. 50 

50. 60 

54.60 

86.80 

96.60 

147. 80 

4 

Percent difference 

-0. 90 

- 0 . 88 

-1. 58 

+0. 84 

-0. 52 

-0. 83 

+ 3. 20 

+ 1. 99 

+0. 30 

+9. 37 


Normalized Modal Mass 

l 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Ml mass configuration (Kgs) 

4. 36 

5. 33 

4. 81 

13.05 

5. 46 

4. 80 

17. 05 

9. 98 

6 . 24 

9.69 

M2 mass configuration (Kgs) 

3.68 

5. 50 

4. 45 

12. 84 

3. 38 

5.66 

20. 17 

6.88 

3.99 

6 . 85 


1 Undamped eigenvalue results. 

2 D2 damping: damped eigenvalue results. 

3 D3 damping: damped engenvalue results. 

4 Based on experimental frequency compared with undamped eigenvalue results. 
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two Booster modes and only one Orbiter mode. The convergence of the 
predicted frequencies is evident as we proceed from the lower to the 
higher number of modes used in the analysis. 

Table VIII(b) shows the natural frequencies and modal masses 
obtained for the free-free component synthesis for the Ml mass con- 
figuration. The highest frequency which could be predicted was for the 
sixth mode with poor results thereafter. This is also evident in the 
mode shapes obtained as shown in Figure 14a. Similar results were 
anticipated for the M2 mass configuration, and thus, an analysis for 
this configuration using free-free component modes was not made. 

B. Damping 

The System damping used to predict System responses was de- 
rived from the method outlined in Section VI. In addition to the damping 
values obtained by our method, we also considered damping values which 
were obtained from a complex eigenvalue solution as discussed in 
Section III. D. 4. These values give an indication of the amount of error 
which could be expected if a straightforward synthesis of damping through 
the coupling matrices was done. Table IX gives the System modal damp- 
ing factors obtained by our method and those obtained from the eigen- 
value results. These are presented for the LI link condition only. As 
pointed out in Section V, component damping energy was obtained ex- 
perimentally for both the pin- slip and free -free configurations. Our 
damping prediction method was applied to both sets of curves, and the 
results are presented in Table IX. Related data in terms of damping 
energy per cycle is presented in Figure 15 for the free-free synthesis 
method. Although the predicted results are shown as smooth curves, the 
individual predicted points deviated from the curve by about the same order 
of error as the experimental points did in the original component curves 
(Figures 8 and 12). 

In general, it appears that the System damping factors obtained 
from the free-free damping energy curves agree with the experimental 
factors somewhat better than those obtained from the pin- slip energy 
curves. The average percent difference is of the order of ±15 percent 
for the damping factors synthesized from the free-free curves and is 
slightly higher for the pin- slip synthesis. This range of error is similar 
to that present in the original component data (Figures 8 and 12), and is 
an extreme improvement over those values obtained directly from the 
eigenvalue results, where the difference between predicted damping and 
experimental damping varies somewhere between 100 and 325 percent. 

To obtain the additional modal damping for the various link dampers 
(L2, L3, L4), we simply add the predicted energy dissipated by the link 
damper to the predicted energy dissipated by the remaining structure. 















































TABLE IX. SYSTEM DAMPING 
(For LI Cases only) 


S-F-Ll -D1 -Ml 

Mode Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 ! 

Energy Prediction* 

.0071 

. 0089 

. 0065 

. 0050 

. 0049 

. 0042 



. 0038 

. 0035 

Experimental Results 

. 0070 

. 0074 

. 0047 

. 0046 

. 0039 

. 0037 

— 

--- 

. 0028 

. 0032 


S-F-Ll -D2-M1 

1 

Mode Number J 

l 

2 

3 

4 

5 

6 

7 

8 

9 

10 j 

Energy Prediction* 

. 0143 

. 0128 

. 0103 

. 0066 

. 0075 

. 0066 



. 0057 

I 

. 0049 ] 

Energy Prediction^ 

. 0138 

. 0100 

. 0088 

. 0061 

. 0060 

. 0053 





. 0046 

. 0040 1 

Eigenvalue Results 

. 0243 

. 0135 

. 0147 

. 0277 

. 0102 

. 0099 

. 0082 

. 0133 

. 0069 

.0103 { 

Experimental Results 

. 0143 

. 0102 

. 0079 

. 0071 

. 0067 

. 0053 




. 004 3 

. 0049 ! 


S-F-L1-D3-M1 

Mode Number i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Energy Prediction* 

• .0296 

. 0237 

. 0196 

. 0130 

. 0133 

. 0112 



. 0087 

. 0078 

Energy Prediction^ 

.0330 . 

. 0226 

. 0190 

.0113 

. 0120 

.0103 



. 0080 

. 0069 

Eigenvalue Results 

. 0562 

. 0278 

. 0311 

. 0607 

. 0189 

. 0166 

. 0133 

. 0256 

.0118 

. 0193 

Experimental Results 

. 0332 

. 0232 

. 0172 

. 0143 

.0115 

. 0105 





. 0080 

. 0086 


S-F-Ll -D2/D3-M1 ' 

— 

Mode Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 ! 10 

Energy Prediction* 

.0169 

.0238 

.0136 

.0102 

. 0088 

. 0071 



. 0067 

.0052 

Energy Prediction^ 

. 0159 

. 0208 

.0113 

. 0099 

. 0074 

. 0060 



. 0059 

. 0047 

Experimental Results 

. 0169 

. 0230 

. 0106 

. 0106 

. 0073 

. 0065 





. 0057 

. 0054 


S-F-Ll -D2-M2 

Mode Number 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


Energy Prediction* 

. 0124 

.0091 

. 0080 

. 0058 

. 0054 

. 0053 


.0041 

. 0044 


Experimental Results 

.0137 

. 0106 

.0071 

. 0071 

. 0050 

. 0056 



. 0067 

- 0056 




1 

2 


Using pin-slip modes 
Using free-free modes 
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The predicted link damping energy is given by D oJ g = TT C0^ el uu, where the 
link damping coefficients C are given in Table III and 0 rel are theo- 
retical mode shape values. Thus, the damping factor for both the struc- 
ture and link damper is the sum of the individual damping factors as 
given by Equation ( 1 1 6 ) . This also provides a further cross-check of 
data by comparing the sum of measured damping factors for the link 
acting alone and oil acting alone, to measured values obtained when both 
act simultaneously. Table X gives the damping factors for the various 
structure damping configurations in combination with the various link 
dampers. 

C. Forced Response 

Forced response measurements were made on the Booster at 
node 10 and on the Orbiter at node 4 for a constant input force level at 
node 10 in the lateral (z) direction. Figures 16, 17, and 18 show the 
experimental and predicted forced responses in the form of response in 
DB versus frequency, similar to the procedures outlined for the compo- 
nent response. They are based on the use of Equation (62), with 
System natural modes obtained from pin- slip synthesis and modal damp- 
ing given by Tables IX and X. Pin- slip, rather than free-free synthesis 
eigenvectors were used in the expansion of forced response, since less 
pin- slip synthesis System eigenvectors are required for sufficient 
accuracy above the sixth mode, as has been indicated in Figure 14. 

The measurements were made on the S-F-L1-D2-M1 and 
S-F-L<1-D3-M1 configurations to determine the accuracy of the mathe- 
matical model and damping prediction for two different structural damp- 
ing conditions that did not include link damping. These are Figures 16 
and 17, respectively. In addition. Figure 18 shows the results obtained 
whereby we added the additional L3 link damping to the S-F-L.1 -D2-M1 
configuration. It should be noted in this figure that the predicted re- 
sponses seem to be consistently higher than the measured responses. 

We suspect that this was a result of a change in the calibration of the 
force link used to drive the structure, but have not verified this. A 
comparison of the predicted responses and the measured responses from 
the above figures indicates that the agreement is very good to excellent 
over the entire frequency range and in particular at the peaks, where 
damping is the dominant factor. 

It must be emphasized that in all forced response predictions, the 
predicted damping factors were used, not the experimental ones. Assum- 
ing that the mathematical model exactly duplicated the dynamic modal 
mass and stiffness behavior of the Shuttle model, then the only difference 
between predicted and experimental responses would come as a result of 
the differences in predicted and experimental damping. Our mathematical 
model, as is the case with all mathematical models, does not exactly 
duplicate the dynamic behavior of the Shuttle model, although judging by 
the frequencies and mode shapes obtained, it is a fairly good one. 



TABLE X(a). SYSTEM DAMPING 
(L2 Case Only) 
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Based on component free-free damping energy curves 
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Based on component free-free damping energy curves 
Obtained directly from decay records 






TABLE X(c). SYSTEM DAMPING 
(L4 Case Only) 
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Based on component free-free damping energy curves 
Obtained directly from decay records 

Obtained by adding the experimental £ due to link damping alone to the experimental £ due to System damping alone 
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Figure 16. System Response for S-F-L1-D2-M1 Configuration 
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Nevertheless, the differences we obtain between predicted and experi- 
mental responses arise now due to the differences in predicted and 
experimental damping, plus the differences in the mathematical and 
real-life Shuttle models. Considering these two factors in combination, 
it is believed that the System forced response results obtained in this 
study indicate the feasibility of substructure testing with subsequent 
synthesis of the component results. 

D. Verification of Damping Synthesis Hypothesis 

It has been pointed out in Section VII, that a necessary part of our 
damping synthesis method includes the hypothesis that experimentally- 
determined curves of the form given in Figures 8 and 12 can be used to 
obtain appropriate damping energies for the components at kinetic energies 
that correspond to points intermediate to the original components. We 
now consider a verification of this hypothesis through the use of System 
damping data. 

Let us consider for the moment, that our primary intent in this 
study was to obtain data that would verify the above hypothesis. We 
start with the curves as they appear in Figures 8 and 12, the points 
for which have been determined by component resonance tests. As has 
been previously mentioned, we do not consider off- resonance testing to 
obtain intermediate values, because of the inherent inaccuracy asso- 
ciated with such measurements. Therefore, in order to obtain values 
at intermediate kinetic energies, and still allow the use of resonance 
testing, we must change the boundary conditions on the respective 
components. At the same time, we must preserve the essential form 
of damping mechanisms in the structure. This is precisely what was 
done in combining the components into the System tests. Therefore, 
we can consider our measurements of System damping, and use these 
data to fill in points on the original component curves. In other words, 
we use the data in the inverse form to which it has been considered 
previously. 

As an example, we consider the free-free component curves in 
Figure 12. We wish to use the System damping results from the various 
configurations to fill in additional points on these curves (assuming the 
data will fall on the respective curves). To do this, we must first note 
that the System damping energies must be properly apportioned to the 
Booster and Orbiter, for each System mode. This is done on the basis 
of the amount of kinetic energy that each component contributes to the 
System mode, which is consistent with the conceptual basis for plotting 
Figures 8 and 12. 

Define average damping energies per cycle as 
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D cs = 4ttT ob (H7a) 

av 

D c0 = 4 ttT 00 G 6 (117b) 

av 

where T o6 and T o0 are the respective Booster and Orbiter kinetic 
energies at a maximum modal amplitude of x 0 , and G s is the measured 
value of System damping in a given mode. Then note that we can write 
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(118a) 

(118b) 
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where D cB1 , 2 an< i D c0l>3 are desired intermediate damping energy 
values corresponding to kinetic energies T oB and T o0 , and both are for 
maximum modal amplitudes of x Q . Further, T o0 corresponds with 
a maximum modal amplitude of x' (we assume that x / <x 0 for a given 
System mode here), and T oS is System kinetic energy. Equations 118 
provide the desired apportioning of damping energy according to the 
respective kinetic energies present in a given system mode, and con- 
tain four unknowns. Two more equations can be formulated from 
Figures 12a and 12b as 


Deal “ 

^cBl 

(119a) 

^cB2 = ^2 

^c02 

(119b) 


where k x and k 2 are the ratios of the curve values at respective 
kinetic energies T qB and T o0 . 


Equations (118) and (119) can be solved to provide 


c 8 av 03 


cBX ~ (T o 8 +KT #0 ) 


C 02 


°av ^'° S 

(T 0 o+ a *2 T oB ) 


(120a) 


(120b) 


Note that if the Booster and Orbiter component curves (Figures 12a, b) 
happen to coincide, then k t = k 3 = 1, and Equations 120 reduce to 
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^cBl 
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CD 

^c02 

= D =o, 


(121a) 

(121b) 


Also note that care must be exercised in the use of Equations (120) so 
that the appropriate amplitudes x 0 and x are assigned in each respec- 
tive System mode. 

System data for several configurations were used in Equation 
(120b) in order to obtain the results shown for the Orbiter in Fig- 
ure 19. Except for one or two isolated points, these data confirm con- 
clusively the original hypothesis, within the same order of error as was 
present throughout the damping measurements. It is particularly note- 
worthy that the data for the D2/D3 configuration also fell on the desired 
curves, since the Booster and Orbiter damping characteristics are more 
significantly different for this configuration. A plot similar to Figure 19 
could also have been produced for the Booster, and comparisons with 
pin- slip configurations could also be employed. However, Figure 19 
should be sufficient to make our point. 
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VIII. FINAL DISCUSSION 


The results of this study have quite dramatically shown the 
successful development of a method for predicting damping in combined 
structural systems with accuracy that is commensurate with that of 
original measurements made on the substructures, for a relatively 
simple parallel- stage Shuttle model. The application of this method 
to a more representative model, and for that matter to a prototype, 
must still be done with caution, and some comments in this regard 
will be given shortly. 

It is first appropriate to discuss the significance of the form of 
the damping energy results in Figures 8 and 12. In essence, each 
curve represents an aggregate or spatial integration of all of the complex 
damping processes which occur in the structure. No assumption is made 
about the particular form of damping present, except that linearity with 
amplitude was found to be present in these experiments. At the same 
time nonlinearity with kinetic energy level exists. We note that the curves 
for pin- slip configurations in Figures 8a, b are somewhat different than 
those for free-free conditions in Figures 12a, b, for corresponding mass 
and damping conditions. In effect, this indicates a significant redistri- 
bution of the damping processes for the different boundary conditions. 

It was determined that this difference results primarily from the pre- 
sence of longitudinal motion in all pin- slip modes, and the form of damp- 
ing chosen, which provided higher effective damping for motion in that 
direction. Of course, this same damping mechanism influenced the 
free-free modes only through rotational displacements, since the damp- 
ing reservoirs were located off the longitudinal axis of the models, and 
only one (Rigid-Body) orthogonal free-free mode includes longitudinal 
displacements. At the same time, it appears that the basic component 
motions (as far as influence damping) in the System modes are more 
like those of the free-free component modes, rather than the pin-slip 
modes, since better damping predictions are provided by free-free mode 
synthesis. A more careful investigation of various modal displacements 
would also shed further light on this behavior. 

The separation of results into groups representing a given aggre- 
gate damping process may cause complications in the analysis of a more 
complex structure. For example, certain kinds of modes may form one 
group (or line on a plot) while other modes may form another line. Some 
evidence of this possibility is indicated in the SAD- 6 data shown in 
Figure 20, which is based on results of experiments given in Reference 8. 
The data are plotted here for fixed amplitude, and as a result, form a 







considerable rearrangement from the form in which they appeared in 
Reference 8. Further, values were obtained by averaging the multiple 
results of modal mass obtained from tests with different suspension 
systems. Also mode numbers are indicated in order of increasing 
natural frequency. We recognize that caution must be exercised in 
interpreting the results as presented here. However, it appears that 
two basic curves are formed at the same modal amplitude, one corre- 
sponding to overall bending modes, and the other to dominantly tank 
cluster motion. Such multiple curves may very well occur in the even 
more complex Shuttle component configurations, and some means of 
handling them will have to be developed. For that matter, some un- 
certainty exists as to what degree grouping of the modes will occur. 

Only analysis and testing of a more representative model, which will 
be accomplished in the extension to the current program, can answer 
these important questions. 

The matter of damping nonlinearity with response amplitude has 
been mentioned in a previous section. The concept of plotting damping 
energy per cycle versus peak kinetic energy is entirely amenable to 
such a complication. So also is the method of damping prediction 
developed herein. However, considerable effort will be required in 
additional testing in order to determine its influence, and it appears 
that at best, much judgement will have to be exercised to include its 
effects in any prototype tests. Such complications may also be the 
subject of extended work. 

It should finally be recognized that the results of this study form 
a far more basic discovery than just prediction of damping in combined 
structures. There is the definite implication that damping energy plots 
can be used to identify similarities in various different complex struc- 
tures. In particular, the slopes, shape, curvature, etc. may eventually 
be cataloged to identify damping processes in many structural applica- 
tions. Thus, each curve is a plot of some now only empirical damping 
law, and the identification of parameters with which these laws can be 
described is an important new development in a very old, yet relatively 
little understood, subject in dynamics. 
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APPENDIX A 


VALUES FOR ELEMENTS OF DISCRETE MODELS 



A- 1 


Orbiter Mass, M t (Kgs) 

Booster Mass, M t (Kgs) 

Mj 

= 1.83 

m 7 

= 1. 83 

Mg 

= 1. 72 

m 8 

= 1. 72 

m 3 

= 1.83 

m 9 

= 1.83 

m 4 

= 1.83 

M 10 

= 1. 83 

m 5 

= 1.40 

M u 

= 1.40 

Mg 

= 1. 83 

M 1S 

= 1. 83 



m 13 

= 1. 83 

Orbiter Moment of Inertia, 
Ii (Kgs-cm 2 ) 

Booster Moment of Inertia 
Ii (Kgs-cm 2 ) 

h 

= 28.78 

I 7 

= 28.78 

h 

= 24.84 

Is 

= 24.84 

h 

= 28.78 


= 28.78 

h 

= 28.78 

*10 

= 28.78 

h 

= 19.20 

Iix 

= 19.20 

h 

= 28.78 

■*•13 

= 28.78 




= 28.78 


Tube thickness, t = 0. 089 cm 
Tube radius, R = 1.18 cm 

Modulus of elasticity, E ■ = 7. 30 X 10 s Newton/cm 2 (10. 6 X 10 6 lb/in 2 ) 
Poisson's ratio, v = 0.325 
\ - 8. 25 cm = 2. 54 cm 


6 3 = 1 1. 43 cm 



A-2 


Distance Between Nodes, i> ifJ (cm) Ml Configuration 

Orbiter Booster 

h* = 24.76 *,. e = 25.40 

>',.3 25.40 i e „ " 25.40 

V = 25 - 40 *•-» * 25 - 40 

= 25.40 Vu * 25.40 

= 24.76 V,» = 25.40 

V.J3 = 24.13 


Distance Between Nodes, X 1)j (cm) — M2 Configuration 

Orbiter Booster 


l ,2 


2,3 


3,4 


■ 4,5 


5,6 


24. 76 

15.24 
30.48 
30. 48 
24. 76 


* 7 , 8 = 25.40 

* 8 , 9 = 10 . 16 

* 9 ,io = 30. 48 

* 10,11 = 35. 56 

*11,12 = l 5 -24 

* 12,13 = 34. 29 



APPENDIX B 


DERIVATION OF ELEMENT STIFFNESS MATRIX 



B - 1 


The derivation of the stiffness matrix for a beam element of length 
& is based on Castigliano 1 s First Theorem. If a structure is subjected 

* * J 1 « ■ 

to a set of external forces P 2 , . . . P n having displacements u 2 , • . . 

u in the direction of these forces, then by the principle of virtual work 
the r-th force may be given by 



‘ (B-l) 


where U is the strain energy of the structure. If we now take a partial 
derivative with respect to displacement u 8 then the stiffness coeffieicnt 
k rfl is given by 


dP r _ 5 s U 
9u s “ du r 3u fl 


(B-2) 


To derive the stiffness matrix for the beam element we first ex- 
press the strain energy in terms of nodal displacement Uj, u 3 , . . . u g as 
shown in Figure B-l 


u s 



Figure B-l. Nodal Displacements of Beam Element I, J 

Following the sign convention shown in Figure B-l, we may write the 
following shear, moment, and slope relationships: 


dV 
d x 

dM 


0 


(B-3) 


dx 


V 


(B-4) 



B-2 


where 


d0 M 
dx = E=9 



-0 


V 

+ nRtG 


(B-5) 

(B-6) 


V = net shearing force on a cross-section 
M = net moment on a cross-section 
0 = measure of net rotation 

E = modulus of elasticity 

G = shear modulus 

=9 = moment of inertia of tube cross-section 

R = mean radius of tube 

t = thickness of tube wall 

u z = displacement in z-direction 

In the axial direction, the displacement u x within the element may 
be written 



These stiffness coefficients are given in the element stiffness matrix. 
Equation (B-24). 


(B- 10) 



B-3 


In the lateral direction we assume a displacement u z in the form 

u z = cl 3 + a 4 x + & 6 x s + a s x 3 (B-ll) 


The coefficients 0- 3 , , a g , and a, g may be found from the condition that 

at x = 0 


u z =u g ; du 7 /dx = -u 3 + 


V(0) 


3 TTRtG 


and that at x = 


u z = u 6 > 




du z /dx = -u, + nKtc . 


(B- 12) 


(B- 13) 


After some manipulation and using Equations (B-4) and (B-5) we find that 
u z can be written as 


where 


u z = u 3 - 


u 


fau.-u,) + u e )} 

1 f J 1 9 J 

+ jr4; — |*3(u s - u 2 ) + l Ui (2A lfJ u 3 +B 10 u 6 )l x 2 
-75“^ [ 2 ( u 5 ~u 2 ) + *> iti (u 3 + u s )l x 3 (B- 14) 


A i»j - 1 + 


3E=9 


B 


i,j 


1 - 


"i., = 1 + 


rrR tG 
6E J 

£? .TTRtG 
i 

12E 

4 l^Rtc 


The strain energy due to bending is given by 


U v 


Ec9 


AJ . . K 

j m 


dx 


(B- 15) 
(B - 1 6) 
(B-17) 


(B- 18) 


while the strain energy due to shear is given by 



B-4 


U s = •§ J XY 2 dx (B- 19) 

0 

where 

X = TTRtG 

Y = slope due to shear (V/X) 

Using Equations (B-4), (B-5), and (B-6) we get 

d 3 u 

V = -E3 £ (B-20) 

dx 

which we substitute into Equation (B-19) to get 

^1,1 V 2 

18E 2 3 s P /d 3 u.\ 

U, = — j— J d* (B-21) 

0 


The total strain energy due to bending and shear is given by 



(B-22) 


From Equation (B-14)we have u z which we use in obtaining U bs . 
stiffness coefficients for bending and shear are now 

ft, r = 2,3,5,6 

rs Bu r 9u a s = 2, 3, 5, 6 

The required stiffness matrix may now be written as 

l>rJ = 



0 

0 - 

ANi,, 

3 

0 

0 

0 

12 

7-5“ 

" i> j 

6 

0 

12 

6 

0 

6 

} 

3+N!,j 

0 

6 

-^i»j 

3-Ni,j 

ANi,,, 

3 

0 

0 

ANi,, 

3 

0 

0 

0 

12 

&t,l 

6 

0 

12 

&t>l 

6 

&ir) 

0 

6 

&Uj 

3-Nt,j 

0 

6 

•^i»j 

3+Ni 
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(B-23) 


(B-24) 



